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Introduc  tion 

Spherical  harmonics  are  closely  associated  with 
developments  of  gravitational  or  magnetic  fields. 
Their  main  importance  in  geodesy  has  its  roots  in 
the  (basis)  property  that  the  external  gravitational 
potential  of  the  earth  may  be  approximated  by  suit¬ 
able  linear  combinations  of  outer  (spherical)  har¬ 
monics  . 

Seen  from  potential  theory  the  standard  represen¬ 
tation  of  a  potential  V  defined  outside  a  sphere 
about  the  origin  with  radius  R  by  a  series  of  outer 
harmonics 

ri 

V(x)  =  ?  (~)n+  I  P*  (cos$)  {A  cos(mX)  ♦  B  sin  (mX)}  (0.1) 

n=o  r  *=o  na  nm  nn 

with 


A  ,  B  :  potential  coefficients  of  degree  n  and 
nm  run 

order  m 


A 


no 


i 

/  V(R,«,X)P*  (cos*)  sin*d*dX 
*  no 

o 


(0.2) 


na 


;  V(»,»,X)  <;$$>  si* 


is  well  -  considered  as  a  series  expansion  by  mul¬ 
tipoles  in  the  center  of  the  sphere,  i.  e.  the  field 
of  a  homogeneous  and  isotropic  distribution  may  be 
adequately  approximated  by  developments  of  spherical 
harmonics.  On  the  other  hand,  the  series  of  spheri¬ 
cal  harmonics  is  less  qualified  for  representing 
the  field  of  a  source  distribution.  Because  of  the 
slow  convergence  of  the  expansion  comparatively  la¬ 
borious  work  must  be  done  to  approximate  local (den¬ 
sity)  anomalies  by  a  partial  sum  of  the  type 


N  -  i  n 

V„(x)  *  I  (-)n+1  Z  P*  (cos*)  iknm  cos(«X) 
N  n=o  r  ®so  nn  nm 


B  _  sin(sX)}. 


(0.3) 


The  price  to  be  paid  is  the  choice  of  a  relatively 
high  degree  N.  This  is  why  in  the  last  decades  many 
efforts  have  been  done  to  reduce  algorithmic  diffi- 


culties  and  computer  time  in  generating  spherical 
harmonics . 

But  not  only  the  slow  convergence  for  modelling 
local  anomalies  proves  a  difficulty  in  using  spheri¬ 
cal  harmonic  expansions.  In  practical  applications, 
almost  all  activities  to  compute  spherical  harmonics 
numerically  are  restricted  to  the  generation  of  one 
and  only  one  system  of  representation  (related  to 
spherical  coordinates),  viz.  the  system 


((i)"1  Y*  .) 
r  J 


n — 0,1,... 

2n+ 1 


(0.4) 


where  Y*  is  defined  as  follows 
n  i  J 


• 

U) 

- 

p* 

(  cos© ) 

n ,  1 

no 

• 

n ,  2 

(O 

= 

P* 

nl 

( cos© ) 

cosX 

n,3 

• 

■ 

(O 

P*  . 
nl 

( cos© ) 

• 

sinX 

• 

■  * 

n ,  2n 

<€> 

. 

P* 

nn 

• 

( cos©  ) 

cos  ( : 

Yn,2n+1(5) 


P*  ( cos©  )  sin  ( nX ) 
nn 


(0.5) 


x  =  r  C,  C 
r  =  |x|  , 


(sin©  cosX,  sin©  sinX,  cos©) 
0  <  X  <  2« .  0  <  ©  <  « 
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Indeed,  as  it  stands,  the  system  (0.5)  has  many  at¬ 
tractive  features.  Perhaps  the  most  important  are 
its  simple  structure  and  recursive  computability. 

A  number  of  different  equation  sets  has  been 
described  in  the  literature  for  suitably  generating 
the  system 


f  y*  i 

n , j  n=0 , 1 ,  . .  . 

j=l i . • • , 2n+ 1 

(e.  g.  Cl.  Muller  (1966),  M.  Gerstl  (1978),  C.  Rizos 
(1979),  0.  Colombo  (1982),  R.  H.  Rapp  (1982).  A  com¬ 
parison  of  methods  has  been  given  by  C.  C.  Tscher- 
ning,  R.  H.  Rapp,  C.  Goad  (1983)). 

In  doing  so  special  care  is  always  needed  to 

ensure  computational  stability  in  the  aeneration  of 
the  associated  Legendre  functions.  It  is  a  critical 

point  to  any  calculation. 

Furthermore,  it  should  be  mentioned  that  a  po¬ 
tential  V  is  related  to  polar  coordinates  when 
choosing  an  approximation  of  type  (0.3).  Thus  dif¬ 
ferentiations  with  respect  to  cartesian  coordinates 

cannot  be  performed  in  a  straiqhtforward  manner.  In 
order  to  determine  the  gravitational  part  of  the 
force  acting  on  an  artificial  satellite  close  to 
earth,  for  example,  we  need  essentially  the  gradient 
of  the  external  gravitational  potential  of  the  earth 
Its  representation  must  be  free  of  singularities 
caused  by  a  special  choice  of  a  coordinate  system. 
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This  is  done  in  canonical  manner  by  use  of  spheri¬ 
cal  harmonic  representations  with  respect  to  carte¬ 
sian  coordinates.  In  addition,  cartesian  coordinates 
seem  to  be  more  adequate  if  the  attempt  is  made  to 
approximate  the  gravitational  part  of  the  earth's 
gravity  potential  by  means  of  outer  harmonics  (mul¬ 
tipoles)  using  a  non-spherical  earth's  model  (e.  g. 
ellipsoid,  spheroid,  telluroid)  as  proposed  by  the 
author  ( 198  3  )  • 

Thus  the  question  arises  of  developing  appro¬ 
priate  algorithms  for  numerical  computation  of 
spherical  harmonics  related  to  cartesian  coordina¬ 
tes  thereby  avoiding  any  computational  problem  of 
accuracy  and  stability. 

The  purpose  of  this  report  is  to  demonstrate 
both  that  spherical  harmonics  can  be  evaluated  exact¬ 
ly  using  exclusively  integer  operations  and  that 
the  procedure  of  expanding  the  external  gravitational 
potential  of  the  earth  into  spherical  harmonics  can 
be  performed  adequately  by  use  of  outer  harmonics  in 
terms  of  cartesian  coordinates. 

The  report  is  meant  to  be  a  proposal  for  exact 

computation  of  spherical  harmonics;  it  is  not  un¬ 
derstood  to  be  the  final  word  in  speed  and  economy 
in  computations  of  e.  g.  high  degree  harmonics. 
Nevertheless,  at  least  in  the  opinion  of  the  au¬ 
thor,  it  is  a  remarkable  discovery  that,  to  any 


r 


prescribed  degree  n,  orthonormal  systems  of  sphen 
cal  harmonics  can  be  deduced  exactly,  i.  e.  exclu- 
sivelv  bv  addition,  subtraction  and  multiplication 


of  integers ,  without  having  any  knowledge  of  spheri 
cal  harmonics  of  orders  different  from  n.  Hence,  as 
long  as  we  (are  able  to)  use  integer  operations,  pr 
blems  of  accuracy  and  stability  do  not  occur;  the 
algorithm  depends  merely  on  available  computer  time 
and  memory. 


Moreover,  our  approach  provides  us  with  exactly 
given  linearly  independent  systems  of  homogeneous 
harmonic  polynomials  (in  terms  of  cartesian  coordi¬ 
nates)  which  can  be  used  for  developments  of  the 
earth's  gravitational  potential  in  series  of  outer 
harmonics  (adapted  to  a  non-spherical  earth's  mo¬ 
del). 

In  detail  we  are  concerned  with  the  following 
considerations : 

Exact  Computation  of  Spherical  Harmonics: 


Starting  point  is  the  class  P  of  homogeneous  har¬ 
monic  polynomials  of  degree  n.  Each  H  e  P  may  be 

n  n 

represented  in  the  form 


iMr 
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(a:  multiindex  of  degree  n).  It  is  obvious  that  the 
set  of  monomials 


forms  a  basis  for  the  space  P  .  The  number  of  such 
monomials  is  precisely  equal  to 

M  =  M(n)  =  (n^2).  (0.8) 

fi 

H  denotes  the  class  of  polynomials  in  P  that  are 
n  n 

harmonic 

H  =  (H  e  P  |  4  H  ( x )  =  0}  .  (0.9) 

n  n  n  x  n 

The  dimension  of  H n  is  2n+l.  In  other  words,  there 

are  2n+l  linearly  independent  functions  H  ,  in  H  : 

n,  j  n 

H  .(x)  =  Z  Cj  x“  .  (0.10) 

n’J  [ a ] =n  “ 

( j=l, . . . , 2n+ 1 ) 

The  main  result  of  the  paper  is  that  the  coeffi¬ 
cients  CJ  can  be  obtained  as  integers  by  integer 
operations  and  that  the  solution  process  of  determi¬ 
ning  the  matrix  c  consisting  of  the  column  vectors 

(C*)  ,  ...  ,  (C2n+1) 
ct  a 

can  be  made  surprisingly  simple  and  reasonably  effi¬ 
cient  . 


In  view  of  the  homogeneity  the  functions  H 


may  be  rewritten  as  follows 


n.  j 


H  .  (x)  = 
n ,  j 


i  n 


[a]=n 


o 


(x  =  |x | 5  ,  |c|  =  1) 


The  functions 


H  .  (C)  =  l  ci  £*  , 


[  a  ]  =n 


(0.11) 

(  j  =  1 , . . . , 2n+l ) 


form  a  maximal  linearly  independent  system  in  the 
space  of  all  (surface)  spherical  harmonics  of  de¬ 
gree  n.  Consequently,  the  system 


( 


2n+l  Hn,j  (x)  *j=l . 2n+ 1  (|x!>0>  (0.12) 


represents  a  maximal  linearly  independent  system  of 
outer  (spherical)  harmonics  of  order  n.  Obviously, 
the  system  (0.12)  does  not  consist  of  polynomials 
but  it  is  a  set  of  algebraic  functions  which  are 
homogeneous  of  degree  -  (n+l). 


Corresponding  to  the  maximal  linearly  indepen¬ 
dent  system  (0.11)  of  surface  spherical  harmonics  of 
order  n  there  exists  a  system 


Approximation  by  aeries  expansion: 


Any  function  V  satisfying  the  following  properties 

(i)  V  is  continuous  in  |x|  R  and  twice 

continuously  differentiable  in  |x|  >  R 

(ii)  V  is  harmonic  in  |x|  >  R 

(iii)  V  is  regular  at  infinity 

may  be  represented  in  the  form 

Zn+1  _  11 

•  &  iii  h  h|/  *;>>  Wi®-  (o.  16) 

(x  a  rC  ,  r  a  | X  I  ,  Id  *  1) 

More  explicitly  this  reads:  given  an  error  bound 
t>  0,  then  there  exists  an  integer  N  =  N  (£)  such 
that 

N  2n+1 

l«  <">  -  A  *  “C"n)  S«.J(n)  ^  S».3<e)l  -  £  <0*17) 

holds  for  all  points  x  e  R-*  with  |x|  ^  R  >  R. 


From  the  addition  theorem  (decomposition  formu¬ 
la)  it  follows  that 

2n+ 1 

v^TTT  P*  (tn)  =  4z.  Y*  U)  Y*  .  (n)  (0.18) 

no  J  =  1  n ,  j  n,j 

2n+  1 

=  .1.  S*  .  (C)  S*  .  (n) 

J=1  n,j  n , j 


for  any  two  unit  vectors  E  ,  n.  Therefore,  for  every 
index  N  and  all  x  c  RJ  with  |x|  _>  Rq  >  R  our  appro¬ 
ximation  to  V 


N  2n*1 


N(x)  =  Z  I  (-pr)n+1  (t-  /  V(Rn)  S*  .(n)  dw}  S*  . 

N  n=o  j=1  f x )  ht  j  i^  n,j  n,j 


(t)  (0.19) 


coincides  with  the  (standard)  approximation 


N  2n*1 


.<«>  ■  &  jfi  (i7T)"*1  hl{,  «*">  *“>  *:,j  <«•  <°-2°> 


Furthermore,  each  element  Y*  .  is  expressible  in 

n,  j  K 

terms  of  the  system  {S*  . ,  and  vice 

n  ,  j  j  =  l . 2n+ 1 

versa.  In  other  words,  even  if  we  base  our  approxi¬ 
mation  process  on  the  conventional  system 


(Y*  . }  „  ,  we  are  able  to  do  this  by  means 

j=l, • • • , 2n+ 1 

of  the  functions  S*  .,  i.  e.  suitable  linear  combi- 

n ,  j 

nations  of  the  system  (0.11). 


12 


The  paper  ends  with  some  reflections  about  the 
role  of  outer  harmonics  (multipoles)  in  orthogonal 
expansions  using  a  non-spherical  earth's  model.  The 
mathematical  procedure  and  the  theoretical  back¬ 
ground  developed  by  the  author  (1983)  are  briefly 
recapitulated.  Many  advantages  of  the  classical 
(strictly)  spherical  approach,  of  course,  cannot  be 
maintained,  when  outer  harmonics  will  be  used  as 
trial  functions  in  approximations  of  the  external 
gravitational  potential  adapted  to  a  non-spherical 
earth's  model  (e.  g.  ellipsoid,  geoid,  telluroid, 
spheroid,  and  (at  least  in  principle)  real  earth). 
Nevertheless ,  numerical  examples  make  us  hope  to 
implement  an  outer  harmonic  (multipole)  approach  to 
the  external  gravitational  potential  of  the  earth 
related  to  a  non-spherical  earth's  surface. 

Seen  from  mathematical  point  of  view,  the  sample 
examples  open  at  the  same  time  new  numerical  per¬ 
spectives  in  solving  exterior  Dirichlet's  boundary- 
value  problems  by  orthogonal  expansion  using  outer 
harmonics . 
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1.  Notations 


3 

R  denotes  three  dimensional  (real)  Eucilidean 
space.  We  write  x,y,...  to  represent  the  elements  of 

3 

R  .  In  components  we  have 


x  =  (x1,x2,x^)T  ,  y  =  (yjtyg.Yj)7  •  (1.1) 

Scalar  product  and  norm  are  defined,  as  usual 


x  •  y  = 

xtyt  +  x 

2y2  + 

x3y  3 

(1.2) 

2 

x  = 

XX  = 

*1  + 

2 

X2  + 

2 

X3 

(1.3) 

M  ■ 

S| 

ii 

/x1  + 

~ 2 - 

X2  + 

—? 

x3 

(1.4) 

Let  a  = 

(al>  *2’ 

a  ) 

be  a 

triple  of 

non-negati- 

ve  integers 

°V  a2  ’ 

a ^ .  We  set 

a !  = 

V  “2! 

°3 ! 

(1.5) 

Ca]  = 

“l  +  °2 

+  “3 

• 

(1.6) 

We  say  a  = 

(al’  “2’ 

xT 

“3) 

is  a  i 

multiindex 

of  degree  n 

if  [a]  =  n,  n:  non-negative  integer. 


As  usual,  we  set 


al  “2  “3 

=  X1  x2  x3 


-  14  - 


(-i-)*1  (-3-)“2  (-9-)°3 

»x/  (  3x2;  lax3; 


(1.7) 


al  “2  °3 

®x  i  ax2  axj 


An  easy  calculation  gives 


(Xj+Xg+X.)"  = 


t  a]  =  n 


n !  a 

— r  x 

a ! 


(1.8) 


(x  y.+x0y„+x-y-)“  =  £  xay01 

11  22  33  [a]=n  “• 


Furthermore ,  we  have 


<»,>**  .  ( 


0  for  a  ^  8  and  [a]  =  [B] 


al  for  a  =  B  . 


(1.9) 


Each  x  e  R  ,  x  =  (x^Xg.x^)  with  |x|  +  0,  ad¬ 
mits  a  representation  of  the  form 


x  =  r5  ,  r  =  |x|  ,  5  =  Uj.Cg.C,) 


)T 

y  ' 


(1.10) 


where  C  e  RJ ,  |C|  =  1,  is  the  uniquely  determined 

directional  (unit)  vector  of  x.  The  unit  sphere  in 
R  will  be  called  0.  As  it  is  well-known,  the  total 
surface  ||ft  ||  of  fl  1st  equal  to  4«: 


|  0  ||  =  /  du>  =  4  * 

n 


(du>:  surface-element), 


12  3 

Let  e  ,c  ,eJ  form  the  (canonical)  orthonormal 


basis  in  R' 


e1  .  0  .  c2  =  ? 


e  =  0 

1 


Then  we  may  represent  the  points  £eft  in  polar  coor¬ 
dinates 

£  =  t  +  i 1-t2  (cosX  e1  +  sinX  e2)  (1.11) 


-  1  <  t  <  +  1  ,  0  <  X  <  2i  ,  t  =  cos  %  , 


(#:  polar  distance,  X:  geocentric  longitude) 


T  T 

£  a  (sin#  cosX  ,  sin#  sinX  ,  cost)  . 


In  terms  of  polar  coordinates  the  Laplace-operator 
A  in 


A  _  „  /  0  v  2  /  0  \  2  /  9  , 2 

A  =  V  •  7  =  (•=— — )  +  (-3 )  +  ( ■  ■ ■  ■  ) 

x  xx  3x  ^  8x2  9x3 

is  represented  as  follows 


,  3.2  23  1  .* 

‘SF*  '  7  IF  *  T  *f  ■ 

r 


(1.12) 


(1.13) 


where  A*  denotes  the  (Laplace-)  Beltrami  operator 
of  the  unit  sphere  ft 


=  d-t2)^)2  -  2t  Ft  +  ^2  (af>2  ‘  <l-“> 


\*  •.*  '  V*  *  •  “  •  »  “  •  '  *  ’  .  *  «  *  «.  *  V  '  ■»  '  »  •  .  •’ .  *.  •*.  «'  •*  *«,*  ’.**. 


If  such  an  operator  is  applied  on  a  homogeneous  po 


lynomial  of  the  same  degree  n 


L  (x)  = 


n 


D0  * 


B 


[  B  3  =n 

we  obtain  as  result  a  real  number: 


3  a2  3  a3  0 

[H  (7  )]L(x)=  Z  Z  C  D  LCl-x1)^  **)(—- J)3 
n  x  n  r  *|  r0T  a  p  a„  1  _  L  a,  i 

W=.  CS]=«  ^ 


Z  C  0  a!  . 

Wa  a 
=n 


Clearly,  we  find 


CV7x)]  Ln(«>  -  [L„(7x>]  »„<*>' 


CH  (V  >]  H  (x)  >  0 

n  x  n  — 


(2. 


(2. 


(2. 


This  gives  rise  to  introduce  an  inner  product 

(•  ,  *)p  on  the  space  P  by  letting 
n 

(H  L)p  =  [H  (V  )]  L  (x)  (2. 

n  n  r  n  x  n 

n 

for  H  and  L  homogeneous  polynomials  of  degree  n. 
n  n 

Clearly,  we  have 


(i) 

(c  H  ,  L  ) 0 
n  *  n  P 

n 

“ 

c  (Hn  •  Ln>P 

n 

for  c 

(ii) 

<VKn-Ln>P 

n 

= 

< Hn ' Ln > P  *  <K 
n 

n’Ln)P 

i 

f or  H  , L  , K  c 
n '  n  ’  n 

P 

n 

(iii) 

<VLn>P 

= 

(L„<H„»p 

(iv)  (H  ,H  )p 
n  nr 


>  0  for  Hn  4  0  . 


The  space  P equipped  with  the  inner  product  ( •  , » )  ^ 
is  a  finite-dimensional  Hilbert  space.  The  set  of  1 
monomials 

{/|7  I  [<*]=«}  (2.8 

forms  a  complete  and  closed  system  of  orthonormal 

elements  in  P  . 

n 


For  each  H  c  P  we  have 
n  n 


H  (x) 
n 


1  aT  [H„<7V>]  y°  x“ 

[  a  3  =  n  a!  n  y 

‘W1  nT  r  5  TT  ^ 


(2.9 


c„„„y)]  iajfl 


Ca]=n 

\ n 

•y> 

n ! 


-4  (x  V  )n  H  (v)  . 


In  other  words 


«„<*>  '  (1^-  '  '  (2',0) 

n 

Pn  equipped  with  the  inner  product  (*,*)p  is  a 

M(n)-dimensional  Hilbert  space  with  the  repro¬ 
ducing  kernel 

K(n;x,y)  =  *  (2.11) 

i.  e.:  (i)  For  every  fixed  y,  the  function 
K( n ; • ,y )  belongs  to  P  . 

(ii)  For  any  e  Pr  and  any  point  x  the 
reproducing  property 


H  (x) 
n 

=  ( K ( n ; x  ,y  )  ,  Hn(y))p 

n 

is  valid. 

K(n;x , y ) 

is  the  only  reproducing  kernel  in  P  . 

Now,  let 

{H  .}  . 
n , j  j= 1 , 

<Ln,j'j=l . N  be  tw° 

or thonorma 1 

systems  in 

P  : 
n 

(Hn  i  • 

n  »  J 

Vk>p  - 

n 

V 

(2.12) 

(L  , 

n  ,  J 

Ln,k)P 

n 

V  • 

(2.12) 


Then,  for  j  =  we  have 

M 

H  T  L  .  (H  .  ,  L  .  )  n 

n,j  .  «  n,k  n,j  ’  n,k  P 

k=  i  n 


"n,  j 


M 

l  Hn  w  .  Hn  w)p 

i_a  n  ,k  Yi ,  j  n  ,  K  r 

k=  l  n 


(2. 13) 


Therefore  it  follows  that 

M  MM 


l  H„  .(x)  Hn  .(y)  =  l  l  L^k(x)  Hr  ,(y)  (H 


j=1 


j=1  k=1 


Ln,k)P 


(2.14) 


M 

I 


L  .(x)  l  . 
n»J  ".J 


M  M 

H(y)  =  1  l  l 

j=1  k=1 


n  i(x)  Hn  k(*>  (L„  i 
n  i  j  n  ^  t  j 


Hn,k}P 


Consequently , 


j=l 


H  .  (x) 
n,  J 


H 


n  i  J 


j  =  i 


n,  J 


.  ( x  )  L 


n,  J 


(y) 


(2.15) 


Hence,  in  particular  for  the  orthonormal  system  of 
monomials  (2.8),  we  obtain 


M 


be  an  orthonormal  system  in 


Let  {Hn  <}i  1 

n,j  j  —  ’  i  •  •  •  i 

P  .  Then 
n 


/  x  n  M 

J  4 


3 

Given  M  points  e  R  and  M  values 

e  R,  we  will  be  able  to  solve  the  inter¬ 
polation  problem 

M 

I  bk  Hn)k(xj)  =  .  j  =  1 »  •  •  •  i  H  ( 2. 16 ) 

k=  1 

if  and  only  if 


det  (H 


n  ,k 


( x  . ) ) 


j=l, 

k=l. 


(2.17) 


A  system  of  M  points  Xj,...,xM  is  called  fundamental 
system  relative  to  Pn  if  the  matrix 


‘V* 

k=  1 , . 


,M 


(  2 . l8  ) 


is  of  maximal  rank  M. 


Let  us  prove  the  existence  of  a  fundamental  sy¬ 
stem  relative  to  F  .  As  orthonormal  system  (H 
the  functions  Hn  M  are  linearly  indepen¬ 

dent.  Hence,  there  exists  a  point  x.  for  which 


Now,  there  must  also  be  a  point  Xg  such  that 


n,l(xl} 

n,2(xl) 

Hn,2<x2* 

4  o 


for  else  we  would  have  a  contradiction  to  the  li¬ 
near  independence  of  H  .  ,  H  In  the  same  way 

Q|  1  n  f  L 

the  existence  of  a  point  x^  can  be  deduced  by  the 
requirement 


that 


Hn, l(xl) 

Hn,l(x2> 

Hn, l(x3J 

«n,2<xl> 

Hn,2(x2) 

Hn,2(x3) 

4  o 

Hn,3(xl) 

Hn,3(x2) 

Hn,  3(x3) 

lly,  we  obtain  a  system  of  points  x^, 

H  , (x,  ) 

e  e  e 

H  ,  (xM ) 

n,  1  1 

• 

e 

e 

n,  1  M 

• 

e 

e 

4  o 

Hn,M(xl) 

e  e  e 

Hn,M(xM) 

i.  e.  a  fundamental  system  relative  to  P  . 

n 


To  every  H  e  P  ,  there  exist  real  numbers 
n  n 


such  that 
M 

H  =  l  b.  H 

n  L  lr 


k  n  ,k 


(2.19) 


Now,  assumed  is  a  fundamental  system  rela¬ 

tive  to  P  ,  the  linear  equations 

M 

Hn,j(xk)  =  bk  ’  k=1 . M  •  (2.20) 

are  uniquely  solvable  in  the  unknowns  a4t...,au. 

1  M 

Thus  we  obtain 


M  M 

l  l  H  (x.)  H 

k=l  j=l  ”  ^ » J  ^  ,  j 


(2.21) 


Let  ^  be  an  orthonormal  system  in 

^n*  Aasume  that  M  is  a  fundamental 

system  relative  to  P  .  Then,  each  H  e  P  is 

n  n  n 

uniquely  representable  in  the  form 


1  a.  K  ( n  ;  x  .  ,  x  ) 
j  =  l  J  J 


i‘- 


3*  Homogeneous  Harmonic  Polynomials 


Let  Hn  c  Pn  be  the  class  of  all  polynomials  in 

P  that  are  harmonic: 
n 

H  =  (H  e  P  I  A  H  (x)  =  0}  .  (3.1 

n  n  n  1  x  n 

For  n  <  2,  all  homogeneous  polynomials  are  har¬ 
monic. 


For  n  >  2,  let  H  „  be  a  homogeneous  polynomial 
■“  n-  a 

of  degree  n-2f  i.e.H  _eP  Then,  for  each 

n  “  a  n*  6 

homogeneous  harmonic  polynomial  Kr,  we  have 


<|x|2  Hn-2<x)  •  ltn(x)) 


(3.2 


[Hn-2  <V3  Ax  Kn(x) 


0 


This  means  |x|2  H  „(x)  is  orthogonal  to  K  (x)  in 

ii"6  n 

the  sense  of  the  inner  product  (*,*)p  . 

n 

Conversely,  suppose  that  Hn  is  orthogonal  to  all 
elements  L  of  the  form 


then  it  follows  that 


0 


(3.3) 
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(  |x  I  2 

H  0(x)  ,  H  (x 

n-2  ’  n 

[H  0 

(  7  )  ]  A  n  ( x  ) 

n-2 

x  x  n 

(Hn-2 

all  H  _  e  P  This  is  true  only  if  AH  =  0, 

n- &  n—  &  n 

e.  is  a  homogeneous  harmonic  polynomial. 


Pn  ,  n  _>  2 ,  is  the  orthogonal  direct  sum  of 

and  Ha  where 
n  ’ 


is  the  space  of  all  members  L  with  L  (x)  = 

n  n 

|x|2  Ln_2(x),  L>n_ 2  e  Pr_2,  i.  e.:  each  homo¬ 
geneous  polynomial  of  degree  n  can  be  decom¬ 
posed  uniquely  in  the  form 

Hn(x)  =  Kr(x)  +  lxl2  Hn-2^x^’ 

where  is  a  homogeneous  harmonic  polynomial  of 
degree  n  and  H  „  is  a  homogeneous  polynomial  of 

IT™  M 

degree  n-2. 


Denote  by  p  and  px  the  projection  operators  in 
onto  W  and  Ha,  respectively.  Then, 


26 


In  other  words, 

K„<x> 

|x|2  Hn_2<„) 

Clearly , 

2 

P  =  P 


p  Hn(x) 

(3.5) 

Px  Hn(x)  * 


(3.6) 


For  all  H  ,  L  e  P 
n '  n  n 


(pH„  , 

ii 

B 

•  PLn 

)p 

n 

(3.7) 

Moreover , 

we  have 

ii 

c 

X 

a 

PKn 

B 

II 

e 

(3.8) 

The  dimension  N(n) 

of  H 
n 

follows  easily 

from  the 

fact  that 

N(n) 

=  dim 

H 

n 

(3.9) 

=  dim 

Pn  - 

dim 

Hi 

n 

=  dim 

Pn  - 

dim 

e 

CM 

1 

*LC 

Explicitly 

,  this  reads 

dim  H  = 
n 


( n+2 ) 
v  2  ’ 


(3.10) 


2n  +  1 


4.  Addition  Theorea 


Our  aim  now  ia  to  give  the  explicit  representa¬ 
tion  of  the  orthogonal  projection  pHn  of  a  given 
homogeneous  polynomial  Hn>  For  that  purpose  we  need 
some  preliminaries. 


A  result  due  to  E.  W.  Hobson  (1955)  yields  for 
i  =  1,2,3 


/  8  .n  1 

V  FI 


(-I)"  (2")i  1 

(n!)  2°  M2"-1 


[n/2] 

[  l 


s=0 


(-D* 


ni(2n— 2s) ! 
(2n)l(n-s)!s! 


where  we  have  set 


[n/2]  =  j  (n  -  j  Ll  -  (-1)"]) 

i  .  e .  : 


n 

2 

for 

n  even 

n- 1 

f  M 

n  o  dd  • 

2 

i  ur 

In  other  words,  we  find 


(e\  )"  Jj 
x  ]7y 


(_Dn  ^2n^! 
(ni)  2n 


[n/2] 

C  I 

s=0 


(-1)S 


n!(2n-2s)! 

(2n)!(n-s)!s! 


|x  |2sAS](e^  *) 
(i  *  1,2,3) 


Since  the  differential  operator  A  is  invariant  with 
respect  to  rotations,  it  is  easy  to  see  that 


<*v 


F 


(-Dn 


(2n)» 
(n!)  2" 


Co/2] 

[  l  (-DS 

s=0 


n ! (2n— 2s ) ! 
(2n)!(n-s)!s! 


|x|2sAS]  (yx)f 


3 

is  valid  for  every  y  e  R  .  Now,  as  we  have  seen  in 
chapter  2,  each  e  P^  may  be  represented  in  the 
form 

H  (x)  =  l  c  (x  x)n, 

j  =  l  J  3 

where  c.,  are  suitable  coefficients  and 

J 

{x.}  is  a  fundamental  system  relative  to 

P  . 

n 


Consequently,  we  have 


Let  be  a  homogeneous  polynomial  of  degree  n. 


Then 


for  x  c  R3,  | x |  4  0 


[H  (7  )] 
n  x '  |  x 


(4.  1 


Cn/Z] 


/  (2n)!  1  r  r  ,  n!(2n-2s)!  ,  , 

<‘,>  77^7  7^ 1  jo  w>  Tmrsns,  w 


x|ZS  AS]  H  (x) 
n 
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Prom  the  considerations  given  in  chapter  3  it 
follows  that 

[Hn(V]  “|7[  =  tKn(V]  ITT  +  [Hn-2(7x)Ux  f^T  *  (4*2> 


Thus,  in  connection  with 


x  TxJ 


A  K  (x)  =  0 

x  n 


we  obtain  for  |x  |  4  0 


CHn(Vx)]  ih  =  [Kn(Vx)]  7^ 


4  o 


(4.3) 


X  £  , 


(4.4) 


.  (.1)“  -iiSiJ - JL _  k  («)  . 

(nl )  2n  |  x|  2n+l  " 


Therefore,  by  comparison  of  (4.1)  and  (4.4),  we  get 


Let  H  be  a  homogeneous  polynomial  of  degree  n. 
n 


Then 


PH„(x) 


trV21/  .xs  n  !  (  2n-2s  )  ! 
l  (-D  ^  i 


I ( 2n-2s  )  !  2s  s  „  f  . 

n  7  M  n-s )  !  s  J  n 


(4.5) 


Observing 


A„(xy)n  =  n(n-l)  | y | 2  (xy)n“2  ,  y  c  ,  (4.6) 


we  obtain,  in  particular 


(4.7) 


Cn/2  ] 

I 

3=0 


/  i\a  (2n-2s) ! (n! )2 

(n-2s) ! (n-s) is ! (2n) ! 


|y|2s(xy)n-2s 


Hence,  after  an  easy  calculation,  we  find  by  using 
polar  coordinates 

x  =  |x|  C 

y  =  |y|  n 

the  equation 

P  (-^f^t  (4.8) 


e  E  Q 

n  e  n 


(2n+l)  2n-n! 


Cn/2] 


( 2n-2s ) ! 


ilLAu  (ixi  iyi)n  i 

1;'  s=0  2n(n-2s)!(n-s)!s! 


(Cn) 


n-2s 


2n+l 


Cn/2] 


1-3’ ... • 2n+ 1 


I x i n i y i n  i  (- i >3  -  (2n-ay).i 


s=0  2n(n-2s) ! (n-a) !s! 


Un) 


n-2s 


Let  {H  . }  . 


be  a  maximal  orthonormal 


n,  j  j= 1 ,  • • • , 2n  + 1 

system  in  H  and  { L_  u  ^  be  a  maximal 

n  n  f  j  j  * 1 1 ■ • • i \ «n+ 1 / 


orthonormal  system  in  Hi.  Then  the  union  of  both 

n 


systems 


(H 


=  2n+ 1 V  ^  Ln , j  ^  j  =  1 , . . . , M- ( 2n+ 1 )  ^  ^ 


SA 


( xy ) 1 

n! 


(4.10) 


2n+l  M-(2n+l) 

l  H  (x)  H  (y)  +  l  L  (x)  L  (y) 

j ^  J  n t  J  j  ^  Hj  j  J 


for  any  two  elements  x,  y  e  R  . 

Furthermore,  in  view  of  the  definition  of  the  pro¬ 
jection  operator  p,  we  get 

2n+l  M-(2n+l) 

p  (  l  H_  Jx)  Hn  i<y>  +  l  Ln  ,(x)  L  (y)) 
n,j  n,j  j-j  n,j  n,j 


2n+l 


M— ( 2n+ 1 ) 


*  P  H»,JU)  Hn,J(y>>+  p  <  ^  Ln,j',t)  Ln,j,y)) 


2n+l 


j=l 


Hn  Hn  * 

n»J  J 


(4.11) 


On  the  other  hand,  as  we  have  shown  above 

in 


P  ( 


(x«y  ) 
n ! 


) 


(2n*1)  2n  n!  i  in  i  in  ,  „\ 
(2n*l) !  K  (_1) 


s  (2n— 2s ) ! 


s=0  2  (n-2s)!(n-s)!s! 


(ft) 


(4.12) 


n-2s 
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By  comparison  of  (4.11)  and  (4.12)  this  yields  the 

addition  theorem  of  homogeneous  harmonic  polynomials 


Let  {H  .  }  .  .  _  .  be  an  orthonormal  system 

n , j  2n+ 1 

in  H  with  respect  to  (•,•)«.  Then,  for 
n  « _ 


x,y  e  R  ,  we  have 


2n+ 1 


(T2 1*1"  l>'npn  <f)  . 

(x  =  |x  |  €  ,  y  =  |y  I  n  ;  £,  n  e  n) 


where  we  have  used  the  abbreviation 


P  (t)  =  Cni23  (-1)3  — - (  -n~2— -  t 

n  a=0  2n( n-2s ) ! ( n-s ) ! s ! 


n-2s 


(t  e  [-1,1]) 


The  function  Pr  ,  n=0,l,...,  given  by 

P  (t)  =  P„A(t)  =  ~  <-!>“  (J)(2n-2s)t2n-2s 


2n  s  =  0 


(t  e  [-1,1]) 


is  the  well-known  Legendre  function. 


P  is  uniquely  determined  by  the  following  pro- 
n 


perties : 


5.  (Surface)  Spherical  Harmonics 


Let  Sn  denote  the  space  of  all  surface  spherical 
harmonics  of  order  n,  i.  e.  the  set  of  all  restric¬ 
tions 


Hn|n  (5.1) 

of  the  members  H  e  H  to  the  unit  sphere  0.  As  H 

n  n  K  n 

is  assumed  to  be  homogeneous,  the  restriction  H  I  ft 

n 

of  H  to  ft  is  given  by 
n 

Hnu)  =  H«(T^T)  ’  ^  ^  o  ,  e  e  ft.  (5.2) 

The  restriction  map 


is  an  isomorphism  of  H  onto  S  i.  e.  (i)  each 

n  n 

element  S  e  S  is  associated  to  one  and  only  one 
n  n  J 

element  H  e  H  (ii)  if  the  elements  S  ,  T  e  S 
n  n  n  ’  n  n 

correspond  to  the  elements  H  ,  K  e  H  .  then  the  li- 

n  ’  n  n  ’ 

near  combination  a  S  +  b  T  e  S  corresponds  to  the 

n  n  n 

element  a  H  +  b  K  e  H  (iii)  different  elements  of 
n  n  n 

S  correspond  to  different  elements  of  H  . 
n  n 

Let  us  consider  S  as  subspace  of  the  space 
2  ^ 

L  (ft)  of  square  -  integrable  functions  on  ft  (equip¬ 
ped  with  the  inner  product  (•,•)  _  ).  Then  there 

L2(  ft) 

are  defined,  for  elements  H  c  H  ,  K  e  H  ,  the  fol- 

n  n  ’  n  n  ’ 


-  36  - 


lowing  two  inner  products 


«,  Kn)  2  =  TZ  /«  U)  K  (  £)  da,  (<•) 

n  n  t2(n)  **  0  ”  " 


(5.3) 


(H  K  )p  =  [H  (7  )]  K  (x). 
n  n  r  n  x  n 


We  now  prove  a  theorem  which  is  of  basic  impor¬ 
tance  for  our  task  of  exact  computation  of  spherical 
harmonics : 


For  H  e  H  K  e  H  , 

n  n  ’  m  m  * 


-  6™,fss§rrr  ‘W3  *»<*>• 


(5.4) 


6  is  the  Kronecker-symbol . 
nm 


Proof.  By  virtue  of  Green's  formula  it  follows  that 
Kn(x)  “  h  £  {-J^J^TKn(y)-Kn(y)  d«(y)(5.4) 

•s  fl 

for  all  x  e  with  |x|  <  lr  where  denotes  the 
derivative  in  the  direction  of  the  outer  normal. 


Therefore  we  find 


CH  <7  )]  K  (x) 
u>  x  n 


(5.5) 


=  h  £  {tH.(V]  7737]-  aT  Kn(y)-I(n(’')  W  [Hs(V3  P7J}  d“(*) 


'*  •  ’  »  *  •'*  .**  /■  •»*'  «*''•  »**  .N  /'  .v  /»  >'r  A  . ‘h  *.  1  _  > 


For  x  4  y  we  have 


tHm<7x)]'P^yJ 


1  ,  -m  (a.)  I  Hm  <X-y) 

=  ,1?  [x-y|  ^n+1 


This  is  equivalent  to 

r  n  1  (2m)  1  Hm  (y-*> 

CHm<7x)]  ‘FyT  -  m!  f  Ix-yl3"*1 


(5.6) 


(5.7) 


Inserting  (5-7)  into  (5.5)  gives 

[H  ( V  ) ]  K  (x) 
m  x  n 


(5-8; 


C  2m)  j  1  f  I  18  ^  ^  §-  *  (y)_K  (y)  JiL 

’  I.-,!2**1 3"»  "  " 


M2”1 


}  du»  (y) 


It  is  easy  to  see  that  for  m  4  n 

«.<V]  Kn(x)lx=0  *  0  <5'9 

and  for  m  =  n 

«„<7X»  K„(x,lx=0  -  tHm<7x)]  K»<x)  =  “WVp,  •  <5'10 

Therefore  we  obtain 


}  dw(y)  (5.11) 


4 f  l  ‘rjss T  S-  7S5-  fTSTT 

0  |y|  y  y  |y| 


2mm! 


TiSTTI  (Hm’Kn>P 
n 


for  m  4  n 


for  m  =  n 


Since  the  normal  derivative  of  K  (x)  and  H  (x) 

n  m 

equal  to 

37  Kn  (rC)  lr=l  *  n  KnU) 

and 

T5?  Hn, 


respectively,  it  follows  that 


1 

■57 


y) 


=  17  /  {n  H  (5)  K  (5)+(m+l)H  (C)K  (5))  du>  (O 

q  m  n  m  n 

=  £L±T7  *  1  /  H  (C)  K  (C)  dui  (£) 

*»  *  m  n 


Thus,  by  combination  of  (5.11)  and  (5.13), 
finally  obtain  the  desired  result. 


are 


(5.12) 


(5.13) 


we 


For  H  e  H  f 

n  n  ’ 

K 

m 

<H  ,K  >  2 
m  n  L2(fl) 

6 

nm 

a 

n 

o »n  is  given  by 

( 2n+l ) ! 

_  i 

"  2n  n! 

—  X 

[H  (V  )]  K  (x)  . 
n  x  n 


=  l*3*...*<2n+1)  * 


To  any  orthonormal  system  in  Hn  with  respect  to 


(-,*) 


n 


^n,j^j  =  l,...  ,2n+l 


(5.14) 


there  corresponds  an  orthonormal  system  in  1.^(0) 

with  respect  to  (  •  ,  •  )  0 

L2(C 1) 

IS*  l 

n.  J  j  =  l, . . . , 2n+ 1 
given  by 


Let  {S 
in  S 

n 

5  ,  n 


*  } 

nj  j= 1 , . . . , 2n+l 
with  respect  to 

e  0  , 


be  an  orthonormal  system 
(*,*)«  .  Then,  for 

l2(0) 


2n+ 1 

I 


j  =  l 


Sn  i(n) 


( 2n+l )  P  (Cn). 
n 


We  discuss  the  most  frequently  used  special  case  of 
the  addition  theorem. 

Let 


c  = 

(  ]fl- t2  cos  X 

,  Tl-t2  sinX 

,  t)T 

n  = 

(  /  1-s2  cos® 

,  / 1-s2  sin® 

,  a)T  . 

Setting  t 

we  obtain 

a  cos®  ,  S 

a  COST  ,  0  <  ® 

,  T  <  «  , 

5  = 

(sin®  cosX 

,  sin®  sinX  , 

T 

cos® ) 

n  = 

(sinx  cos® 

,  sint  sin®  , 

»T 

cost) 

(0  £  X,  ®  <  2t ) 


Thus  it  follows  that 


For  our 


orthonormal  system  in 

we  pick  { Y*  , 

( Q)  11  »  J  J  “  1  » 


S  with 
n 

• • f 2n+ 1 


respect  to 
(cf.  (0.5)) 


^Ti  Pn  0(t) 

-  p;,0^> 

^(awl)fel7tPnl(,)  “»* 

=  P*  ,(t)  cosX 
n,  l 

/2(2“*l>  tSttt  p„i(t)  •*** 

=  P*  .  (t)  sinX 
n,  l 

/2(2n*,)  -WTT  Pm,<t)  c“(ni) 

/2(2“*1>  ISsH  Pn»(t> 


=  P*  (t)  cos(nX) 
n,n 

=  P*  ( t)  sin(nX) 
n,n 


Then,  by  virtue  of  the  addition  theorem,  we  obtain 


2n+ 1 

l  K  ,(5)  y*  .  (n) 


j  =  i 


n,  j 


n  i  J 


=  (2n+l)  P  (cos®)  P  (cost) 

n  n 

+  2(2n+l)  Z  |  P  (cos®)  P  (cost)  cos  m(X-4>) 

m=  x  v  n+m )  :  rui  rm 

On  the  other  hand,  we  have 


P_  (5n) 


P  (sin®  sinT  cos(X-<|/)  +  cos®  cost) 


Thus  it  follows  that 


P  (sind  sinx  cos  ( X—  4* )  +  coat  cost) 
n 


(5.1< 


P  (cost)  P  (cost) 
n  n 

+  2  7 £ll!! It  P_(cos4)  P_ (cost)  cos  m  (X-4) 

nt  1  v  rwm )  2  nm  m 


which  is  the  classical  addition  theorem  (decomposi¬ 
tion  formula)  for  Legendre  functions. 


Conventionally  the  functions  P*  defined  by 

nm 


p;0<^ 


p*  (t) 
nm 


^^Tpn0(t) 


l'2(2-1>  tS^Stt  p™(t) 


(5-1' 


and  the  trigonometric  expressions 


sin  (mX) 
cos  (mX) 


(5-H 


are  computed  by  recursion. 

The  generation  of  sin(mX)  and  cos(mX),  m  >  2, 
can  be  done  through  the  following  recursion  rela¬ 
tionships  : 


sin(mX)  =  2  cosX  sin  ((m-l)X)  -  sin  ((m-2)x) 
cos(mX)  =  2  cosX  cos  ((m-l)X)  -  cos  ((m-2)x) 


These  relationships  are  useful  for  point  calcula¬ 
tions,  but  they  are  inefficient  for  use  if  a  set  of 
points  at  a  uniform  longitude  interval  is  being 
considered . 

The  generation  of  the  functions  P*  is  "criti- 

nm 

cal  to  any  calculation  involving  spherical  harmonic 
expansions.  In  choosing  an  algorithm  one  must  con¬ 
sider  the  speed  and  stability  and  accuracy  of  the 
procedure."  (R.  H.  Rapp  (1982)).  The  essential  point 
is  the  sensitivity  to  small  computational  errors 
when  n  increases.  In  the  past  few  years  a  number  of 
different  ways  have  been  described  in  the  litera¬ 
ture.  A  comparison  of  methods  has  been  given  by 
C.  C.  Tscherning,  R.  H.  Rapp,  C.  Goad  (1983). 


6.  Exact  Coaputation  of  Linearly  Independent 

Syateaa  of  Homogeneous  Haraonic  Polynomials 

Our  purpose  now  is  to  explain  the  theoretical 
background  of  an  alternate  method  for  computation  of 
spherical  harmonics. 

The  concept  is  based  on  the  observation  that  any 
linearly  independent  system 


^n,j^j  =  l,...  ,2n+l 


(6.1 


of  homogeneous  harmonic  polynomials  of  degree  n 


r  i  c1 

[ a  j  =n  a 


(6.2 


Hn,2n+l(x) 


I  C2n+1  x° 


C  a] =n  ° 


can  be  generated  by  exact  computation  of  the  coeffi- 

cients  ,  j  =  l,...,2n+l  (i.e.  exclusively  by  in- 
a  1 ) 

teger  operations).  That  is, the  functions  H  .  are 

n ,  j 

available  globally  in  closed  form  as  linear  combi¬ 
nation  of  monomials. 


1)  We  write  briefly  instead  of  when  confu- 

a  n ;  a 

sion  is  not  likely  to  arise. 


Given  a  homogeneous  polynomial  Hr  of  the  form 


H  (x)  =  I  C  x“  ,  n  >  2  (6.3) 

n  [o]=n  “ 

Assumed  H  is  harmonic,  i.  e., 
n 

A  H  (x)  =  0  ,  x  e  R3,  (6.4) 

x  n 


we  obtain 


A  H  (x) 
x  n 


A  I  C  x 
X  r  -l  « 
[ot]  =  n 


(6.5) 


Z  C  A  (x  )  =  0  . 


[  a]=n 


Thus  it  follows  that 


a„-2  a,  a,  a.  o.-2  cl  a  a  a,-2 

Z  Cjo^o^Dx^  x2  ’VVVl)x1  *2  *3  ^(«3-1)x1  «2  x?  ](6.6) 

[a]=n 


0 


We  discuss  the  terms 


in  more  detail.  Every  term  in  (6.7)  with  index 
T 

a  =  ( a  ^ , a^ , a^ )  satisfying  +  ag  +  a^  =  n  is  a 
homogeneous  polynomial  of  degree  n-2.  Hence,  the 
sum 

A  H  (x)  (6.8) 

x  n 

o.-2  a 

■rf  C.  C“l<a1-1),1  *2" 

LaJ*n 


3  **2(V1)x1  *2  *} 


a„  a, 


S(v1>,i1 


l2°3 
x2  x3 


-t 


is  a  homogeneous  polynomial  of  degree  n~2.  There¬ 
fore  AHn  can  be  represented  in  the  form 

A  K  (x)  =  I  Dr  x8  .  (6.9) 

X  "  [ B]=n-2  8 

The  coefficients  Dfl  are  given  by 

P 


(a]=n 


C  m. 
a  Ba 


(6.10) 


where  m_  is  defined  as  follows 
oa 


m 


Ba 


\  «!<«!-»> 

for 

B-a  =  (-2,0,0) 

1 

CM 

a 

CM 

a 

f  or 

B-a  =  (0 , -2 ,0 ) 

.3u3-d 

for 

B-a  =  (0,0, -2) 

(6.11) 


0 


otherwise 


H  is  assumed  to  be  harmonic,  i.  e.  A  H  (x)  =  0 
identically  for  all  x  c  R  .  But  this  means  all 
numbers  Dg  are  equal  to  0. 

Hence ,  we  have 


1  c„  =  0 
[ a ]=n  ° 


for  all  8  with  [8]  =  n-2. 


(6.12) 


Now,  (6.12)  is  a  linear  system  of  (”)  equations 

in  (n*2)  unknowns  C  ,  [a]  =  n. 

A  a 


The  matrix 


m  =  (mfla)  (6.13) 

has  (”)  rows  and  (n*2)  columns,  m  can  be  parti- 


,  ,nN  ,  ,n+2 

has  (^)  rows  and  (  2 

tioned  as  follows 


}  (£> 


(6.14) 


(n2)  ( n+22)-(n2 ) 

=  2n+ 1 


where 

1  = 

( 1B6 

and 

r  = 

(rB6 

matrix . 

For 

the 

set  of 

(rB6)  i.  a  (”)  by  <"f)  -  (”) 


For  the  set  of  multiindices  of  degree  n  we  in- 


troduce  a  binary  relation  between  elements 


<*'  =  ,  “2  ’  apT 

a"  =  (a”  ,  0"  ,  o»)T 

designated  by  "  >  "  and  defined  as  follows: 

o'  >  o'»  (6.15) 

if  and  only  if  one  of  the  following  relations  is 
satisfied 

o'  >  o'* 

or 


=  0L\ 


or 


a 


2 


-  ai 


*2 


a2 


“3 


>  o! 


The  binary  relation  "  >  "  implies  an  ordering 
for  the  multiindices  a,  [a]  =  n  according  to  the 


(n,  0 

,  0) 

(n-  1 , 

1. 

0) 

(n- 1 , 

o, 

1) 

(n-2, 

2, 

0) 

(n-2 , 

1, 

1) 

(n-2, 

0, 

2) 

1  )  1 
2  1 


(0,  n,  0) 

« 

(0,  0,  n) 


/  n+2  v 
(  2  )_n 


(n+2) 
V  2  1 


In  the  same  way,  the  set  of  multiindices  8, 

[&]  =  n-2,  may  be  ordered  by  increasing  integers  i, 
1  i  <_  (g).  Hence,  in  canonical  manner,  each  pair 
(B,«)  with  [8]  =  n-2,  [a]  =  n  corresponds  uniquely 
to  a  pair  (i,j),  1  _<  i  <_  (g)  *  1  <  j  ("g2). 

In  this  notation  the  matrix 


m  =  (m^)  ,  [8]  =  n-2,  [a]  =  n 


can  be  rewritten  in  the  ordered  form 


m  =  (m^)  ,  1  <  i  <  (g)  ,  1  <  j  <  ( ”g2 ) 


Analogous ly 


1  =  (lB>f)  ,  [8]  =  n-2  ,  C-rD  =  n-2 


becomes 


1  = 

<1iJ)  •  'i 

i  <  (“) 

•  1  < 

1  <  (S 

From  (6. 

11)  it  can  be 

deduced 

that 

=  0  for  i  > 

j  ,  i 

* 

s 

e 

e 

CM 

II 

C  CM 

xiJ 

4  0  for  i  = 

J  .  i 

—  1  ,  .  .  .  , 

(")  . 

But  this 

shows  that  1 

is  non-s 

ingular , 

hence 

matrix  m  is  of  maximal  rank: 

rk(m)s(g).  (6.l6 

Therefore  we  are  able  to  find  (n*2)  -  (£) .  i.  e. 
2n+l  linearly  independent  solution  vectors 

(A*) . <A2n+1)  ,  [a]  =  n, 

ct  a 

of  the  homogeneous  linear  system  (6.12). 


According  to  standard  conlusions  in  Linear  Al- 
n+2 

gebra  the  (  g  )  by  2n+l  matrix  a  consisting  of  the 
vectors  (A1)  ,  ...  ,  (A2n+1) 

CX  OL 

«  =  (  (  aM . (A2n+1)H  ("*2)  =  M  (6.17 

o  a  « 


2n+ 1 


may  be  partitioned  in  the  following  form 


•  •  o  ■■ 
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’  u]  }  (n22)-(2n+l) 
a  =  -i  }  2n+ 1 

L 

where  i  is  the  (2n+l)  by  (2n+l) 

Then  the  linear  system 
ma  =0  (6.19) 

can  be  written  as  follows 

lu  =  r  .  (6.20) 

Since  1  is  a  (2n+l)  by  (2n+l)  upper  triangular 
matrix,  the  unknown  matrix  u  can  be  computed  by 
(2n+l)- times  backward  substitution. 

The  elements  of  the  matrix  m  are  exclusively 
integers.  Therfore,  any  solution  A^,  [a]  s  n  of  the 
linear  system  (6.19)  is  a  column  vector  of  rational 
components.  Hence,  there  exists  a  matrix 

c  =  ( (C  * ) , . . . , (C2nn ) )  ,  [a]  =  n, 

the  elements  of  which  are  exclusively  integers  (ob¬ 
serve  that  if  (A  ),  [a]  =  n,  is  a  solution  of 

OL 

(6.19),  then  (C  )  =  k  (A  ),  [a]  =  n,  k  integer,  is 
a  a 

a  solution ,  too ) . 

In  other  words,  the  solution  process  can  be 
performed  strictly  in  the  modulus  of  integers. 


(6 . l8 ) 

unit  matrix. 


Exact  computation  (without  rounding  errors)  is  pos¬ 
sible  in  integer  mode  by  use  of  integer  operations 
(addition,  subtraction,  multiplication  of  inte¬ 
gers)  . 

When  the  matrix  c  has  been  calculated,  the  ho¬ 
mogeneous  harmonic  polynomials  H  .  given  by 

n ,  j 

H  .(x)  =  I  x“  ,  j  =  l,...,2n+l 

n ,  j  r  i  a 

LaJ  =  n 

form  a  (maximal)  linearly  independent  system  in  Ho¬ 
llaing  spherical  coordinates 

x  =  | x |  (sin©cosA  ,  sin©sinX  ,  cos©) 

we  obtain 

i  o.,  <*o  01 1 

E  CJ  (sin  ©cos  AXsin  acos 

Co]=n  “ 

The  functions 

(H  l 

n, j f j»l, . . .  ,2n+l 
given  by 

.  (at+a2)  a  o„  a. 

H  .(£)  =  E  CJ  sin  ©cos  Asin  Acos  J© 

niJ 

aj+Og+a^n  12  3 

T 

(£efl  ,  C  =  (sin©cosA  ,  sin©sinX  ,  cos©)  ) 


Hn,j(x)  =  M 


form  a  linearly  independent  system  of  surface 
spherical  harmonics  of  order  n. 
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L*> 


i 

i 

* 

* . 

.* 

S 

i 

S 


r.' 


.  •' 


f-' 


Computational  Aspects:  In  each  row  the  matrix  m 
contains  three  elements  different  from  zero:  the 
first  of  them  is  a  diagonal  element.  Therefore,  the 
matrix  m  (incl.  its  associated  "index  matrix")  re¬ 
quires  only 

|  n  (n-1) 

locations  of  storage  if  an  optimal  storage  scheme 
is  used.  (Remark  that  complete  storage  needs 

n4  -  2n3  -  n2  -  2n 

locations). 

If  we  look  at  the  special  structure  of  the  so¬ 
lution  matrix 

c  =  ( (C1 ) . (C2n+1)), 

01  Oi 

2 

we  have  at  most  n  -n+1  non- vanishing  elements  in 

each  column  vector  (C^).  Indeed,  the  total  number 

01 

is  less.  This  gives  rise  for  a  further  reduction  in 
storage  requirements. 

Finally,  it  should  be  emphasized  that  exact 
computation,  i.  e.  addition,  subtraction,  multipli¬ 
cation  in  integer  mode  must  be  performed  strictly 
in  the  available  range  of  the  integer  constants 
(determined  by  the  architecture  of  the  computer). 


■,*  ».*  */  *, 

■A  ^1  rt 


->  -> 


•*. “V  - « .**. 
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That  means  all  data  generated  in  the  course  of  the 
calculations  should  be  as  small  as  possible.  An 
efficient  tool  is  decomposition  into  a  product  of 
prime  numbers. 

Reaark:  The  loss  of  precision  using  real  computa¬ 
tion  has  not  been  investigated  by  the  author. 

Exaaples :  Let  us  demonstrate  the  technique  of  cal¬ 
culating  the  matrix  c  for  two  examples. 


Exaaple  1:  n  =  3 


Then  an  elementary  calculation  yields 


( n+2 ) 

K  2  1 

<") 

'2 

(n+2)  _  ( n }  _ 

v  2  '  v  2; 


10 

3 

7 


Every  homogeneous  harmonic  polymial  e  may 
be  represented  in  the  form: 


H3*X*  =  C300  X1  +  C210  xlx2 


+  C201  xlx3 


+  C120  X1X2  +  C"4  x'Xr 


111  X1X2X3 


+  C102  xlx3 


+  C030  x2  +  C021  X2X3 


+  C. 


012  X2X3 


+  C003  x3 


(X  =  (Xj.Xg.X.)  ) 
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I 


has  to  fulfill  the  differential  equation 
Ax  H3(x)  =  0  , 

i .  e . 


6  C300  X1 

+ 

2  C 

210  X2 

+  2 

C20 1  x3 

+ 

2  C 120  X1 

+ 

6  C, 

030  x2 

+  2 

C02 1  x3 

+ 

2  C 102  X1 

+ 

2  C 

< 

012  x2 

+  6 

C003  x3 

= 

0 

Since 

Ax  «3(x)  « 

0 

identically 

3 

for  all  x  e  we  get 

(n)  = 

3  equations 

6  C300  + 

2 

C  120 

+  2 

C  102 

=  0 

2  C2 10  + 

6 

C030 

+  2 

C0 12 

»  0 

2  C20 1  + 

2 

C021 

+  6 

C003 

=  0  . 

r.  Using  the  introduced  order  for  the  coefficients  C 

r  Qt 

[a]  =  3 ,  the  equation  me  =  0  reads  in  matrix  nota- 
P  tion 


■» 


r 

k 


^  1  **  »-*  f-~-  »r*  *t" 


a. 
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60012020000 
020 :ooo6ooo 
002:0000206 


c_ 

»  ' 

0 

300 

C210 

0 

C201 

0 

C120 

0 

cm 

0 

C102 

0 

C030 

0 

C021 

0 

C012 

0 

C003 

0 

.  . 

where  we  have  marked  the  partitioning  of  the  matrix 
m  and  the  vector  (CQ)  by  dashed  lines. 


If  we  choose 

C 120  =  -1  *  Clll  =  **•  =  C003  =  ° 
the  linear  system  is  uniquely  solved  by  the  vector 

(-j  ,  0  ,  0  :  -1  ,  0  ,  0  ,  0  ,  0  ,  0  ,  0)T 

Multiplying  this  vector  by  3  all  components  become 
integers 

(c*)  =  (1  ,  0  ,  0  :  -3  ,  0  ,  0  ,  0  ,  0  ,  0  ,  o) 


In  the  same  way  we  generate  a  set  of  7  linearly  in¬ 
dependent  solutions  of  the  above  system  the  compo¬ 
nents  of  which  are  all  integers,  viz. 
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«:> 


=  (o  ,  o  ,  o  :  o  ,  -i  ,  o  ,  o  ,  o  ,  o  ,  o) 


(c3) 


=  (i  ,  o  ,  o  :  o  ,  o  f  -3  ,  o  ,  o  ,  o  ,  o) 


<c«> 


=  (0  ,  3  ,  o  :  0  ,  0  ,  0  ,  -1  ,  0  ,  0  ,  0) 


<C*> 


=  (o  ,  o  ,  i  :  o  ,  o  ,  o  ,  o  ,  -i  ,  o  ,  o) 


«6.> 


=  (o  ,  i  ,  o  :  o  ,  o  ,  o  ,  o  ,  o  ,  -i  ,  o) 


(c7) 


=  (o  ,  o  ,  3  :  o  ,  o  ,  o  ,  o  ,  o  ,  o  ,  -i) 


Thus  the  following  linearly  independent  system 
{H  -i  ^ -s  -i  n  of  homogeneous  harmonic  polynomials 

of  degree  3  is  generated  by  the  following  functions 


H3,1(x)  * 


3  2 

1  •  xJt  -  3  xa  x2 


h3>2(x)  = 


-1 


X1  x2  x3 


H3 . 3 < x  >  * 


X1  ”  3  X1  x3 


H3,4<x)  » 


2  3 

3  xl  x2  “  1  *  x2 


H3,5U)  " 


2  ,  2 
X1  x3  ”  1  *  x2 


H3,6(x)  = 


2  „  2 
X1  x2  -  1  '  X2  X3 


H3,7(x)  = 


2  3 

3  -  1  •  x' 


r„*  Xj  x2  X3 

I  H3,3(«>  -1300 

-3102 


H34(«)  = 


3 

-1 


X1  x2  x3 


2 

0 


1 

3 


0 

0 


H3,5(x)  " 


XXX 

l  2  3 


l 

-l 


2 

0 


0 

2 


1 

1 


<2  )  =  21 


<2>  -  10 


Hence , 


(n+2) 
v  2  ’ 


<2>  =  11 


Every  homogeneous  polynomial  of  degree  5  may  be 
written  in  the  form 


Hj.  ( x ) 


[«]=5 


C  x 
a 


“l  °2  “3 

2  Cn,  «  «  X1  X3 

a1+a2+a3=5  al°‘2°t3  2  3 


Explicitly , 


H5(x)  =  C5(X)  xt 


+  c4l0  X1  x2  +  C401  X1  x3 


+  C320  X1  x2  +  C311  X1  X2  X3  +  C302  X1  x3 
+  C230  X1  x2  +  C211  X1  X2  X3  +  C212  X1  x2  X3 
+  C203  X1  X3  +  Cl4 0  X1  x2  +  C131  X1  X2  X3 
+  C122  X1  x2  x3  +  C113  X1  X2  x3  +  C104  xl  x3 
+  c050  x2  +  C04l  x2  x3  +  C032  x2  x3 


+  C023  x2  x3  +  C0l4  x2  x3  +  C005 


'  — i  ^  ^  ^  »  -  *  —  * 


Hence,  the  equation 


Ax  H5(x)  =  0 


x  e  R' 


is  equivalent  to 

.3  ,  2 


20  C500x1  *  12  C410X1X2  *  12  C4oix1  x3  *  6  C320x1x2  +  6  C311X1X2X 

2 


3 


+  2  C320x1  *  6  C230x1x2  *  2  C221x1  x3  +  12  C140X1X2  +  6  C131X1X2X3 


3 


♦  2  Cj02x1  ♦  2  C210X1X2  *  6  C203X1  *3  *  2  C122x1x2  *  6  C113X1X2X 

2 


6  C502xlX‘  ♦  2  C2J0  t\  ♦  2  C211  x2Xj  ♦  2  x.,x2  ♦  2  x* 


*  2  C122X1X3  *  20  C050  x2  *  12  C041  X2X3  +  6  C032  X2X3  +  2  C023  x3 

*  12  C10Vlx5  +  2  C032  x2  *  6  C023  X2X3  *  12  C014  X2X3  +  20  C005  x3 


=  0  . 


But  this  gives 


20 

C500 

+ 

2 

C  320 

+ 

2 

C302 

12 

c4io 

+ 

6 

C230 

+ 

2 

C2 12 

12 

c4oi 

+ 

2 

C221 

+ 

6 

C203 

6 

C320 

+ 

12 

C 140 

+ 

2 

C  122 

6 

C311 

+ 

6 

C 131 

+ 

6 

C113 

6 

+ 

2 

^  i  o 

+ 

12 

C  A 

0 

0 

0 

0 

0 

0 


+  20  C 


+  2  C 


0 


2  C230 
2  C211  + 
2  C212  + 
2  C203  + 


12  C 


6  C 


2  C 


050 

04l 

032 

023 


o32 


6  C 


+  12  C 


023 

0l4 


+  20C 


005 


0 

0 

0 


These  equations  can  be  rewritten  in  matrix  notation 
as  follows 


20  0  0  2  0  2 
0  12  0  0  0  0 
0  0  12  0  0  0 
0  0  0  6  0  0 
0  0  0  0  6  0 
0  0  0  0  0  6 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 


0000:0 

• 

6020:0 

0206:0 

• 

0  0  0  0  :i2 

• 

0000:0 
0000:0 
2000:0 
0200:0 
0  0  2  0:0 
0002:0 


0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  2  0  0  0  0 
6  0  6  0  0  0 
0  2  0  12  0  0 
0  0  0  0  20  0 
0  0  0  0  0  12 
0  0  0  0  0  0 
0  0  0  0  0  0 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
2  0  0  0 
0  6  0  0 
6  0  12  0 
0  2  0  20 


'500 


'410 


'401 


'320 


'311 


'302 


'221 


'212 


C203 

£*•* 

U140 


'131 


'122 


'113 


'104 


'050 


'041 


'032 


'023 


'014 


|>5J 


=  ( 


where  we  have  again  marked  the  partitioning  of  the 
matrix  and  the  vector  in  the  indicated  way. 


We  choose  C  l4(J  =  -1  ,  CJ31  =  ...  =  CQ05  =  0.  Then 
the  linear  system  is  uniquely  solvable  by  the  the 
vector 


(-  1  ,  o  ,  o  ,  2  ,  o  ,  o  ,  o  ,  o  ,  o  ,  o  :  -i  ,  o  ,  ...  ,  o)T 


Multiplying  this  vector  by  the  factor  5  all  compo¬ 
nents  become  it*  *  egers 

(c1)  =  (-i,  o,  o,  io,  o,  o,  o,  o,  o,  o  :  -5,  0,  ....  o)T 


In  the  same  way  we  generate  a  set  of  11  linearly  in 
dependent  solutions  of  the  above  system  the  compo¬ 
nents  of  which  are  integers  by 

(i)  choosing  the  lower  part  of  the  vector 
identically  0  besides  one  component 

(ii)  solving  the  system  by  backward  substitution 


(iii)  multiplying  every  resulting  vector  by  an 
appropriate  integer. 


According  to  this  procedure  the  following  system  is 
generated 
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7.  Exact  Computation  of  Orthonormal  Syatema  of 
Homogeneous  Harmonic  Polynomials 

Corresponding  to  the  linearly  independent  system 
of  homogeneous  harmonic  polynomials  of  degree  n 

/  u  l 

n ,  j  j= 1 , . . . , 2n+ 1 
an  orthogonal  system 

/  LJ  l 

n, j  j  =  l ,  • • • , 2n+ 1 


can  be  constructed  as  usual  (according  to  the  well- 
known  Gram-Schmidt  process): 


The  functions  H  .  are  computed  recursively.  We 

**  J 

start  from 


H  .  =  H  . 

n, 1  n, 1 


(7.D 


Then  we  set 


H 


n ,  2 


a  ‘  H  .  +  H  . 
2,1  n , 1  n , 1 


(7.2) 


n 


The  coefficient  a„  „  has  to  be  chosen  such  that  H 

2,1  n , 2 

is  orthogonal  to 


(H„,2  •  Hn,l>P  ‘  0 

’  n 


(7.3) 


It  turns  out  that 
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2,1 


(Hn,2  ’  Hn,l)Pn 
(»n,l  ’  »n,l>P 


(7.4) 


It  should  be  noted  that  numerator  and  denominator 
may  be  determined  exactly  (cf.  (2.5)). 


Now,  let 

fi  i  9  4  +  a"  H  _  +  H  „ 

^  i  ^  ^  i  *  ®  i  *  3 »  2  n  ,  2  n  f  3 

The  requirements 


(7.5) 


(H  ,  , 

finfl>P 

n 

=  0 

(H  _  , 
n,3  ’ 

«„,2>P 

n 

=  0 

1  to 

n 

(H  , 
n » 3 

■  fin,l>P 

n 

*3,1  " 

(N 

n ,  1 

■  Vi>p 

n 

n 

( H_  o 

n,  3 

'  fi„,2>P 

n 

*3,2  - 

( N  „ 
n ,  2 

’  ^n ,  2 }  P 

(7.6) 


(7.7) 


Again,  the  coefficients  can  be  deduced  exclusively 
by  integer  operations. 


Analogously  we  get,  in  general, 


H  ,  =  a”  H 
n,k  k,1  n,l 


+  •  •  •  + 


a.n  ,  ,  H  +  H  . 

k,k-l  n,k-l  n,k 

( k  =  2 , . . . , 2n+l ) 


H  =  H 

n ,  1  n ,  1 


where  the  coefficients 


n 


(H  ,  H  )p 
n  ,  K  n  y  s  r 


k ,  s 


(H  ,  H  )„ 
n, s  n , s  F 


(7.8) 


(7.9) 


(7.10) 


are  computable  exactly  by  integer  operations,  i.  e. 

a"  is  known  exactly  as  fraction  of  integers. 

A,  S 


According  to  the  orthogonalization  scheme,  each 

function  H  .  is  a  linear  combination  of  the  func- 
n,  J 

tions 


H  .  .  .  .  .  ,  H  n  .  . 

n,l  *  *  n,2n+i 

The  coefficients  of  this  linear  combination  can  be 
obtained  exactly  as  rational  numbers,  too.  Thus 
there  exists  a  vector 

<cJ> 

such  that 

H  .(x)  =  Z  Ci  x“  ,  j=  1 , . . . , 2n+ 1  .  (7.11) 
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The  vectors  (C^) ,  j= 1 , . . . , 2n+l ,  form  a  matrix  c  ex¬ 
clusively  consisting  of  fractions  of  integers  (as¬ 
sumed  all  numbers  in  the  course  of  computation  have 
been  calculated  in  such  a  way  that  numerator  and 
denominator  are  known  as  integers). 

There  exists  a  sequence  of  homogeneous  harmonic 
polynomials 


{Hn, j}j=lf... ,2n+l 


wi  th 


‘n.j  '  Vl>f_  ■  0  for  •>  +  1  • 


viz . 


H  =  a"  4  H  .  +  . 
n  ,k  k  ,  1  n  ,  1 


.  +  a.  .  ,  H  +  H 

k , k— l  n | K" i  n j K 


(k  *  2 , . . . , 2n+l ) 


H  =  H  . 

n  ,  1  n ,  1 


In  addition,  a  sequence 


{K  .}  i  . 

n,j  j= 1 , . . . , 2n-  1 


(7.12) 


can  be  constructed  so  as  to  have  the  property 


(Kn,j  ’  Kn,l)Pn  =  6jl 


1  for  j  =  1 
0  for  j  ^  1 


The  functions  K 

_  n.J 

each  element  H 

n,  j 


can  be  obtained  by  dividing 
by  its  norm 


n 


H 


n,  J 


H 


n,  J 


The  sequence 

^Kn,j*j=l,... ,2n+l 

forms  an  orthonormal  system  of  homogeneous  harmonic 
polynomials  of  degree  n. 


Provided  the  expression 


P 

n 


has  been  stored  as  radicant  of  an  integer,  the 
functions 


{ K  . }  .  _  , 

n , j  j=l, . . . ,2n+l 

are  available  exactly,  too.  The  exactness  is  of  basi 
importance.  It  implies  the  stability  of  the  solution 
process.  In  fact,  the  method  is  constructed  so  as 
to  have  no  sensitivity  to  computational  errors. 


To  any  orthonormal  system  of 

polynomials  {K  . 

n , J  j— l,..., 


homogeneous  harmonic 
in  (Hn  ,  ( • , • )p  ) 


2n  +  1 


there  corresponds  an  orthonormal  system  of  spherical 

harmonics  <S*  j=1 . 2n+1  in  <«n  ,  <  ’ .  *  >  t2 ( n ,  > 

given  by 


S* 

n 


* 


j 


n 


(7.13) 


According  to  our  orthonormalization  process, 

each  element  S*  .  is  available  as  follows 
n,  j 


S*  ,(5)  =  Z 
n’J  [a]=n 


a 


K  e  0 


where  we  have  set 


Bv 


B 


.  ,  fi  )p 

n'J  n • J  n 


(7.14) 


(7.15) 


Hence,  the  radical  sign  and'  the  division  are  the 
only  sources  for  rounding  errors. 


Since  spherical  harmonics  of  different  order  are 
orthogonal  we  finally  find: 

\  The  system 


.  ,H 

tj’  " 


z 

Qt14Ct2+ot3=n 


V 

“l°2a3 


(a^p  o 
sin  c  ®  cos 


I  2  3 

Xsin  acos  v 


T 

(j  =  l,...,2n+l  ;  5  =  (sin®cosX  ,  sinOsinX  ,  cos®)  ) 


is  orthonormal  in  the  sense 


n,  1 


Sk  2 
k,J  £2(n) 


=  6 


nk 


Summarizing 

Each  element 


our  results  we  obtain: 


S*  .  of  the  the  orthonormal  system 
n ,  J 


{S*  } 

n , j  j= 1 , . . . , 2n+l 

of  surface  spherical  harmonics  of  order  n  may 
be  represented  in  the  form 


S*  (?)  =  z  B-J  c  ,  c  e  n  , 
n’j  [«4.n  “ 


i  (a1^2)  °1  a2  a3 

Z  sin  ©cos  Xsin  Xcos  v 

a1+a2+0t3=n 


(j  =  l,...,2n+l,  £  =  (sinOcosX  ,  sinftsinX  ,  cos©)1)  , 
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where  the  coefficients  are  given  as  follows 

a 


and  the  values 


(fln,J-Hn,J>Pn  *■'“  *1 

can  be  determined  exclusively  by  integer  opera¬ 
tions  . 


Examples:  Again,  we  discuss  the  orders  n  =  3  and 
n  =  5. 


Example  1:  n  =  3 


According  to  the  aforementioned  orthonormaliza 
tion  process  due  to  Gram  -  Schmidt  we  are  able  to 
deduce  from  the  maximal  system  of  linearly  indepen 
dent  homogeneous  harmonic  polynomials 


.7 


an  orthogonal  system 


,7 


The  resulting  functions  are  listed  below: 
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Thus,  the  integers  are  decomposed  into  a  (positive) 
integer  times  a  product  of  prime  numbers  <  3 • 


Consequently,  the  orthonormal  system 

tK3,j}j=l,...,7 

corresponding  to 

>  i  - 

may  be  listed  as  follows: 

Table:  n  =  3 

2  3  xj  x2  x3 

K3,1^X^  = 

100300 
-10  112  0 
131:  normalization-factor 

K3,2*x^  = 

10  0  111 
10  0 

K3,3(x)  = 


100300 
10  0  12  0 
12  0  10  2 


f: 
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Finally,  the  orthonormal  system 
{S*  } 

3, j'j=l,...  ,7 

of  surface  spherical  harmonics  of  degree  n  (with  re 

spect  to  (•,•)  is  given  as  follows 

l2(fi) 

S3.j*^  =  ^“3  K3,j  .  (x  =  |x|  (  ,  ?  t  D) 

wi  th 


=  105  =  1  *3  *5  *7 

Example  2:  n  =  5 

Analogous  calculations  give  the  following  functions 
Table:  n  =  5 


2  ^  ^  X1  X2  X3 

K5,.(*)  * 

-1  0  0  0  5  0  0 

110  13  2  0 

-iooii4o 

1  7  1  1  :  normalization-factor 

K5,2^x* 

1  0  0  0  3  1  1 

-1  0  0  0  1  3  1 


12  10 
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K5 ,3(x) 


K54(x) 


3  5 


X1  x2  x3 


0 

1 

3 

0 

3 

7 


0 

0 

0 

1 

1 

3 


0 

0 

0 

0 

0 

0 


5 

3 

3 

1 

1 


10  0  0 
10  0  0 
110  0 
12  2  0 


0 

2 

0 

4 

2 


0 

0 

2 

0 

2 


3  1  1 
1  3  1 
1  l  3 


K5,5(x)  = 

-1  0  0  0  5  0  0 

-1  1  0  0  3  2  0 

1  2  1  0  3  0  2 

-1  0  0  0  1  4  0 

12  10  12  2 

-1  3  0  0  1  0  4 

7  6  2  0 


-  8l  - 


3  5  xt  x2  x3 


K5,H(x)  = 


-1  0 


4  0 


3  0  12  0 

0  110  4 

3  0  10  2 

3  0  0  0  0 

6  3  1 


The  functions  {S*  }.  now  are  given  as 

j  — 

follows 


s5,j<5)  =  ^5  K5 .  j  ’  (x  =  lxl  5  ^  e  n 


with 


«5  =  10395  =  1.3.5.7.9.11  . 


A  list  of  the  first  orthonormal  systems  of  homo¬ 
geneous  harmonic  polynomials  is  given  in  Appen¬ 
dix  2.  In  the  same  way  as  illustrated  above,  the 

coefficients  and  the  normalization-factor  are 
a 

split  into  an  integar  times  a  product  of  prime 
numbers  <  n. 


o.  Relations  to  the  Standard  Systen  of  Spherical 


Ha  monies 


*Hn,  j  * j  =  l, . . . ,2n+l 


(n  >  2) 


(8.1) 


be  the  sequence  of  linearly  independent  homogeneous 
harmonic  polynomials  of  degree  n  generated  by  exact 
computation  as  described  in  chapter  6: 


H_  .  ( x )  = 


C  a  ] =n 


Cj  x“ 
a 


x  e  R' 


(8.2) 


where  the  column  vectors  (C^)  of  the  matrix 


c  =  ((C*)  .  (C*n+1)) 

a  a 


(8.3) 


consist  exclusively  of  integers.  (Remember  that  the 
matrix  c  has  the  form 


,  p  M-( 2n+l )  , 

q  2n+l  , 


(8.4) 


2n+l 


where  M  is  equal  to  (g  )  and  q  is  a  (2n+l)  by 
(2n+l)  diagonal  matrix  with  integers  different  from 
zero  on  its  diagonal). 

It  is  clear  that  the  system 


-  «  *  *  **>  l.  •*,  '."k  ■' 


given  by 


H 

n, 


m 


Z  cj  5“  ,  5  e  Q 

[«]=n  “ 


Hn,2n+1(5)= 


Z 

[  a  ]  =  n 


C2n+1 

a 


C  c  n 


(8 .( 


forms  a  maximal  linearly  independent  set  of  (sur¬ 
face)  spherical  harmonics  of  order  n. 


On  the  other  hand  we  know  that  the  standard  sy¬ 
stem 


Y»,.<E)  * 

p„0(cose) 

Yn,2<E>  ' 

Pni ( COS« ) 

cosX 

Yn  1<E> 

n ,  3 

• 

• 

P  .(cos©) 
n  1 

• 

• 

sinX 

(8. 

• 

Y„,2„<E>  - 

• 

P  (cos©) 
nn 

cos ( nX ) 

Yn,2„*l(E)  * 

P  ( cos©  ) 
nn 

sin( nX ) 

(Wtj, c2,e3>T  , 

C^sin&cosX  f 

52=si«®sinX 

,  £y=COS©) 

is  a  maximal  linearly  independent  system  of  (sur¬ 
face)  spherical  harmonics  of  order  n. 


Our  aim  now  is  to  point  out  some  relations  between 


these  two  systems  of  linearly  independent  spherical 
harmonics . 


As  maximal  linearly  independent  systems,  both 
are  bases  in  the  space  of  (surface)  spherical 
harmonics  of  order  n.  Thus  each  element  of  the  first 
basis  may  be  represented  by  a  linear  combination  of 
the  second  basis,  and  vice  versa. 


In  particular,  we  have 


*  Mf0  •  5  *  °  <8-8> 

with  suitable  coefficients 


*0 . ®2n 

This  means  that  the  system  (Y  . .  _  .is 

Q  y  J  J  =  1  |  •  •  •  y  «11+  1 

representable  as  sum  of  monomials 


Yn  iU)  =  1  El  * 

n,J  (o]=n  “ 


C  e  0  .  (8.9) 


Therefore,  by  comparison  of  (8.6),  (8.8)  and  (8.9), 
we  obtain 

Y  .(5)  =  Z  Ej 

n  y  j  r  i  a 

w  L  a J=n 


2n  _ 

=  i  gJ  z  c2n+1~M  ea 

y=0  ^  [ a ] =n  “ 


z  zn  gj  c2n+1“^  5“ 

C  a ] =n  u=0  *  “ 


(8.10) 


I  j  c2n+l-n  =  Ej  (8.11) 

_  °u  a  a 

H=0 

(Coi]  =  n  ,  j=l, . . . ,2n+l)  . 

In  vectorial  form  this  yields 

c  g  =  e  ,  (8.12) 


where  c,e,g  are 

given  as 

follows 

rl 

c2 

_2n+ 1 

n,0,0 

n,0,0 

•  •  • 

n  ,0 , 0 

C1 

c2 

r2n+l 

n-1,1,0 

n-1,1,0 

n-1,1,0 

c1 

c2 

r2n+l 

n- 1 ,0 , 1 

• 

n- 1,0,1 

• 

•  •  • 

n- 1,0,1 

C  "  2,n-2,0 

• 

• 

• 

c2 

z,n-2,0 

• 

•  •  • 

• 

• 

„2n+l 
z , n-2 ,0 

• 

• 

• 

C1 

• 

• 

c2 

• 

• 

2n+l 

2 , 0 , n- 2 

2 , 0  ,  n-2 

•  •  • 

2,0, n-2 

(n+1) 

1  2  ' 


( 2n+l ) 


1 , n- 1 , 0 
0  . . 


'l.n-2,1 


iivr 


JWf 


>JfV 


n- 1 , 1 ,0 

n-1,  1,0 

•  •  • 

n- 1,1,0 

1 

n-1 ,0, 1 

E2 

n-1 ,0 , 1 

•  •  • 

-2n+l 
n- 1,0,1 

• 

• 

• 

1 

2 , n-2 , 0 

• 

• 

• 

E2 

2,n-2,0 

•  •  • 

• 

• 

• 

-2x1+1 

2,n-2,0 

• 

• 

1 

2,0 ,n-2 

2 , 0 , n-2 

•  •  • 

• 

• 

• 

_2n+l 
2,0, n-2 

o 

1 

e 

H 

1 , n- 1 , 0 

•  •  • 

1  ,n- 1 ,0 

1 

1 ,  n-2 , 1 

E2 

l,n-2,l 

•  •  • 

_2n+ 1 

E1 ,n-2 , 1 

• 

• 

• 

1 

1,0, n-1 

1,0, n-1 

•  •  • 

• 

• 

_2n+ 1 

E1 , 0 ,n- 1 

K 


According  to  the  ordering  of  the  multiindices  in¬ 
troduced  in  chapter  6  we  are  able  to  rewrite  the 
matrices  c  and  e  as  follows: 


[<; 

c2 

^  ^  •  •  • 

r2n+l 

S 

C2 

*“2  *  '  * 

r2n+l 

C2 

c2 

v>  ^  •  •  • 

r2n+l 

C3 

• 

• 

ci 

CM- ( 3n- 1 ) 

• 

• 

• 

c2 

cM-(3n-l) 

• 

• 

• 

„2n+l 

UM-(3n-l) 

• 

• 

• 

c1 

M-(2n-l) 

• 

• 

c2 

m-(2n+l) 

• 

• 

r2n+l 

LM- ( 2n+ 1 ) 

r  * 

0 . 

.  .  .0 

M- ( 2n ) 

• 

0  . 

• 

• 

• 

C2 

^M- ( 2n- 1 )  *  *  . 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

*  * , 

• 

*  .  • 

'  0 

• 

% 

0 


0 


r2n+l 


From  the  special  structure  of  the  matrix  c  we  de¬ 
duce  that 


S 


2n 


sn+ 1 


n 


g 


0 


i .  e , 


JM-( 2n) 


j  =  1, 


'M-  (  2n ) 


17  J 

M-(n+l) 


j  =  1. 


'M-(n+l ) 


El  n 

M— n 


,n+l 


=  1 


'H-n 


eh 


.2n+l 


j  =  1. 


'M 


,c(2n+l)-nrl 

M-JJ 


(p  =  0  ,  . 


In  order  to  calculate  these  values 


coefficients  E>* 


M- ( 2n ) 


»  •  •  •  »  E 


M 


J  = 


To  this  end  we  recall  to  our  mind 


p  (cose)  =  (sine)m/2  i  (_1)k  (2n"2k)! 
rm 


k=0  2^ !  ( n-k ) !  ( n-m-2k ) 


. . . , 2n+ 1 , 

. . . , 2n+ 1 , 

. . . , 2n+ 1 , 

•  ■ • i 2n+ 1 1 

(8.1 

,  2n) 

we  need  the 
1 , . . . , 2n+ 1 . 

tha  t 

-  (cose)n"m“2k 


-  90  - 


Furthermore,  it  is  well-known  that 

t?  k 

cos(mX)  =  Z  ^cosX)m'2k  (sinX)2k 

K— U 


[—3 
L  2  J 


(8.15) 


,  .  „  -  (-1)  It!  ,  ,  sm-2k-l 

sin(mX)  =  Z  (zk+l)!(m-2k-77T  (c°sX>  (sinX) 


Thus,  an  elementary  computation  gives 


P  (cos®)  cos(mX) 

n* 


[±±]  [i] 

1  *  j  (-1)k+1 _ w!  (2n-2k)! _ 

|x|n  k=0  1=0  2°k ! ( n-k ) ! ( n-m-2k ) ! ( 21 ) ! (m-21 ) ! 


I  1 2k  m-21  21  n-m-2k 

M  *1  x2  x3 


(8.16) 


P  (cos®)  sin(mX) 


c~]  [~] 

J-  z  z  (-i)k+1 


|x|n  k=0  1=0 


m ! ( 2n— 2k ) ! 

2% !  ( n-k ) !  ( n-m-2k  ) !  (  21+1) !  (s-21-1)  1 


I  | 2k  m-21-1  21+1  n-m-2k 
|x|  x1  x  2 


According  to  the  polynomial  theorem  we  obtain 


I  I  2k  ,  2  2  2xk  _  k!  2a 

|x|  =  (x  +x  +x  )  =  *  ZT  x 

L  a  J  =k 


(8.17) 
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P^eos®)  cos(aX) 


[^pl 

=  Z  6 
k=0 


n,« 

k 


1 

(il(k-p)T 


^+2C?  ^+2C?  c;-*-2m  ♦ ... 


(®  ~  0,1,..., a) 


reap. 


(8.19) 


P__(cos6)  sin(aX) 


C^p] 


Z  6«.«  o»  l  _L_  c«-1-2[=p]t2t=p]*2M  +  _ 

_0  k  [l'1]  h=o  P!(k~P'T  ^  c  y 


(a  — 


It  ia  eaay  to  see  that 


-  -n.m  .  y  1  _m-2[™] r2p+2[™]  rn-m-2p 

„!o8k  TC=J  T^TTEhmTT  2t2  S 


rMIl 
L  2  J 
I 

v=0 


z 

k=v 


1 

1 

v! (k-v) ! 


m  jn-2[5]  _2v+2[?]  rn-m-2v 

V]  S  2  *2  2  S 

2J 


(m  *  0,(1), n) 


and 


(8.20) 
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£ 

ftn,m 

k=0 

[5=2] 

L  2  J 

'[5=2] 
L  2  J 

£ 

£ 

v=0 

k=v 

» 

,n,m  1 
k  v! (k-v) 


1  m-l-2[^]e2C^]+2M  cn-m-2M 

c-u) !  i  2  3 


m  r2v-2[^i]  _n-m-2v 

W1  2  '2  2 

L  2  -1 


(m  =  1,(1), n) 


Consequently  we  find 


■  =  0: 


Yn  1U)  = 

n,l 


Pn0(cOs6)  = 


2n 

£ 

p=0 


1 


with 


1 


[§] 

2 

j 

rf1*0  1 

{c(2n*l)-Mr* 

k-Ji 

K~2 

^  (|)!(k-|)! 

m-m 

for  p  =  0, (2) ,n- 

•c^ 

0  otherwise 


1  <  n  <  n^n  even: 


Y  _  (E)  =  P  (cosO  cos(mX) 
n,2m  rm 


2n 

£ 

|j=0 


2n+l-y 


(O 


with 


m 


v  £  8k 
2  kJ^ 


n,m 


<Hg")  !  (k  -  J^S)! 


(c(2n»l)-M(-l 

M— 


t  •«  \  ^ 

for  |j  =  m,(2),n-[l  + - ^ - •] 


0  otherwise 


1  <  ■  <  n;s  odd: 


Y  __(£)  =  P  (cos®)  cos(mA) 
n,  an  nm 


2n 


2m 


I  H„  swi  „<5> 

M=0  ^  n’2n+1_M 


with 


kJi^a 


-n.in _ .  1 _ 

k  (Jd=2=S)  s  (k  -  li=f=!l)  I 


for  ii  =  n+m,  (2)  1  + 


0  otherwise 


1  <  ■  <  n;a  even: 


Y 


n.2m+l 


U)  = 


2 n 
l 


2n+l 


H»,2»*l-„<5) 


P  (cos6)  sin(mX) 


,_(2n+l) 


2m+l 

= 


<T  o 
hh2 


£ 

1 

M-m-n 
K“  2 


e; 


,n,m 


_ 1 _ 

(J£J£),(k-l£*a)i 


for  i<  *  n+m,  (2)  ,2n-[  1 


t-n"*1 

2 


0  otherwise 


1  <  ■  <  n;n  odd: 


Y  _-(C)  =  P  (cos«)  sin(mX) 
n,2m+l  nm 


2n 

I 

M=0 


2n+l 


Hn,2n+1-M(E) 


with 


a n+1 


of" 


rn-m-i 
L  2  J 

Z  8n,m 


(c; 


2n+l-Hj-l 


2 

/  l)n 

for  n  =  m,(2),n-[l  + 


0 


otherwise 


Furthermore  we  get  (cf.  (0.5)) 


2n 


-  '/STT  ^  *1  Hn.2n.lVe> 


1 


-  96  - 


Yn  i(?)  = 

n,  j 


I  (n-C4]) !  2n  . 

2(2n+l)  - f -  Z  gJ  H  _  .  JO 


(j  =  2,(l),2n+l) 


Therefore  the  addition  theorem  may  be  rewritten 


as  follows: 


Let  {H  . .  _  4  be  the  sequence  of  linear- 

ly  independent  homogeneous  harmonic  polynomials 
of  degree  n  generated  by  exact  computation  (cf. 
chapter  6).  Then,  for  any  two  £,n  t  0, 


P„Un> 


2n  2n 


pio  v=0  ***  ^  Hn»2n+1-M(5)  Hn,2n+l-\ 


2n+l  (n-C-i]).'  2n  2n 

+  2  Z  - 7 -  Z  Z  gJ  gJ  H  „  .  (C)  H 

j=2  (n+t4])!  p=0  v=0  ^  V  »«2n+1-M 


n,2n+l-v 


Examples 


We  demonstrate  some  relations  between  the  system 
(8.6)  and  (8.7)  for  the  degrees  n  =  3,5. 


Example  1:  n  =  3 


P^q(cos«) 

= 

^3  2  *1*3  2  *2*3 

= 

Y3.1,E> 

P^1(cos«) 

coaX 

= 

6 

*1*3  2  *1  2  *1*2 

= 

v32(E) 

P^Ccos©) 

ainX 

= 

3 

2  3  r2  3  r3 

52^3  "  ?  5152  “  f  ^2 

= 

Y3.3(E) 

P^2(cos«) 

coa { 2X ) 

= 

15 

«?«3  “  «  C^3 

= 

Y3.4<E> 

P^2(cos©) 

ain(2X) 

= 

30 

w3 

= 

Y3.5<5> 

P^fcoa©) 

coa  ( 3X ) 

= 

15 

^  -  45  *^2 

= 

Y3,6(5) 

P^Ccoa©)  ain(3X) 

s 

W  -  15  C l 

= 

Y3,7<E) 

(£=(£« ,C„i£o)  »  C,=ain©coaX  ,  t_=sindsinX  ,  E_=coa©) 


s  s 
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Moreover,  we  get 


5,.(E> 

= 

^7 

<-* 

H3,5m  ' 

"3.7(C>) 

3,2(E) 

= 

<  1 

H3 , 1 <E )  - 

2H3,3(C) > 

;,3<e» 

= 

r~' l 

(  ? 

"34(C)  - 

3H3,6<  C ) » 

;,4<e» 

= 

rl& 

(  15 

"35(C)) 

5.5<E) 

= 

(-30 

H3i2(5)) 

3,6<e> 

s 

(  15 

"3,(0) 

3,7(E) 

= 

(  15 

"34(C)) 

• 

Example  2:  n  =  5 

Y5.1<5)  ■  4  H5,ll<e)  +  4  H5,9(5)  *  «*  "s,?*5’ 

Y5,2«>  ■  4  H5,5<E)+s8  *4oH5,1(E) 

y5.3<E)  *  4  H5,io(5)  *  4  ‘S,8<5>  *  *5  H5.6<5) 
Y5,4(E>  '  4  H5,9  <E)  +  4  H5,7(E) 

Y5,5(E)  '  4  "5,4  <E>  *  4  ,%.2(E> 

Y5,6<E>  -  4  H5,3  (E>  +  *10  H5,l(t) 

Y5,7(E>  =  *3  H5,8  (E)  *  *5  H5,6(E) 

Y5,8(e>  -  t\  »f7  (C> 
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Y5,9  (E) 

-  4 

H5,2 

W5) 

10 

=  g10 

H5,i 

Y5,.l<5) 

11 
=  S5 

H5,6 

(5) 

U) 

(5) 


We  omit  the  explicit  calculation  of  the  coeff 
cients . 


9*  Representation  of  External  Gravitational 
Potential  by  Means  of  Spherical  Harmonics 
(Exterior  Dirichlet's  Problem) 


A  surface  S  c  RJ  will  be  called  regular,  if  it 
satisfies  the  following  properties: 

3 

(i)  S  divides  three-dimensional  Euclidean  space  R 

uniquely  into  the  bounded  region  E  =  E^  (inner 

space)  and  the  unbounded  region  E  (outer 

3  -  _  e 

space)  defined  by  E  =  R  -E  ,  E  =  E  \j  S .  The 

J  e  e  e  e 

origin  belongs  to  E^  . 

(ii)  S  is  a  closed  and  compact  surface  with  no 
double  points. 


(iii)  S  is  a  C 


-  surface. 


From  the  definition  it  is  clear  that  all  (geo- 
detically  relevant)  earth  models  are  included.  Re¬ 
gular  surfaces  are  for  example  sphere,  ellipsoid, 
spheroid,  telluroid  or  real  (regular)  earth's  sur¬ 
face)  . 

Let  V  be  a  function  satisfying  the  following 
proper ti es : 


(i)  V  is  continuous  in  E  =  E  u  S  and  twice  con- 

e  e 

tinuously  differentiable  in  E 


( || S || :  total  surface  of  S).  According  to  the  expan¬ 
sion  theorem  developed  by  the  author  (1983)  the  po¬ 
tential  V  can  be  represented  by  the  series 


V(x) 


•  2n+l 


(9.5) 


where  the  numbers 


(  V  ,H 


i2(S) 


=  1Jf]f  {  V(y)  Hn,j(y)  du* 


are  the  Fourier  (or  orthogonal)  coefficients  of  V  on 

S  with  respect  to  the  system  {H*  . }  -  . 

n | j  n=o t  *  9 • • • 

j  =  1  ,  e  e  e  ,  2ll+ 1 


More  explicitly  this  reads:  given  an  error  bound 
£  >  0,  then  there  exists  an  integer  N  =  N(£)  such 
that 

I V( x )  -  VN(x)|  <  £  (9.6) 

with 

N  2n+l  . 

Vx)  =  ^  £  {lTTTl/  v(y)  Hn  i(y)  dw}  Hn  i(x)  (9.7) 

iteO  j=l  II  bll  S  n,J  n’J 

holds  for  all  x  e  G  C  E@  and  dist  (G,S)  >  p  >  0.  In 
each  compact  subset  K  c.  E0  with  dist  (K,S)  ^  p  >  0 
the  convergence  is  uniform 

sup  | V(x)  -  V  (x) |  <  £  . 
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The  series  guarantees  ordinary  pointwise  conver¬ 
gence,  in  fact,  in  the  (whole)  outer  space  E^. 

In  the  context  of  boundary- va lue  problems  of  po¬ 
tential  theory 

N  2n+ 1 

Vn<x)  =  Jo  J,  f  1R  l f<y>  d“1  H".J(X) 

may  be  interpreted  as  approximation  to  the  exterior 
Dirichlet's  problem 

AV ( x )  =  0  for  x  e  E 

e 

I  V(x  )  |  =  0  (-ply)  ,  |  7V(x  )  I  =  0 

|X|  |x|2 

( |x I  -  “) 

V(x)  =  fix)  for  x  e  S  (f  e  C(0)(S)) 

by  the  orthogonal  expansion  in  terms  of  outer  har¬ 
monics  H*  .  of  order  <  N. 
n.J  - 

When  our  expansion  theorem  is  formulated  especial¬ 
ly  for  a  sphere  S  about  the  origin  with  radius  R,  we 
are  led  to  the  classical  results  (cf.  Introduction). 


Examples:  In  order  to  get  an  impression  of  the  me¬ 
thod  of  computing  a  potential  V  by  the  the  (gene¬ 
ralized)  Fourier  procedure  proposed  here  we  discuss 
some  sample  examples. 

(i)  Boundaries 

Three  boundaries  will  be  used  in  the  examples, 
their  cross-sections  in  the  (xt  ,  x^)-plane 
are  shown  in  the  figures 


Surface  1 


Surface  2 


Surface  3 


The  first  surface  (Surface  1)  is  a  sphere  gi- 

3  T 

ven  by  x  e  R  ,  x  =  (x^Xg.x^)  with 


Xj  s  *  sin®  cosX 


x2  ~  "5  sin®  sinX 


X  ^  =  ~  cos® 


o  <_  «  _<  k 

0  <  X  <  2* 
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(ii) 


The  second  surface  (Surface  2)  is  an  ellipsoid 
given  by 


X1 

s 

1  sin® 

cos  X 

x„ 

4  sin® 

s  inX 

0  <  « 

<  « 

2 

2 

X, 

_ 

2  cos® 

0  <  X 

<  2 1 

3 

» 

The  third 

surface 

(Surface  3) 

is  given 

by 

X1 

= 

<p 

U 

CM 

sin®  cosX 

X2 

= 

3  r  ( ®  ) 

sin®sinX 

<*> 

v| 

O 

<  1 

x. 

__ 

r  ( ®  ) 

cos® 

0  <  X 

<  2 1 

3 

with 

r(»)  =  (cos  (26)  +  (-ji  -  sin2(2®))1/2}1/2. 


Potentials 


We  want  to  give  expansions  for  the  following 
two  potentials: 


v<1)(x)  >K('' 


3 

|2  /  X2  ,  f  l2 

1  cos  ( - *)  +  e 

|x|2 


sin  (- 
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[(TT2  -  5)2*<t%  -  -  3)2J1/2 

1*1  1*1  1*1 


(x  =  (x^x^xj) 

The  potentials  ,  V^2^  are  the  Kelvin 

transforms  of  the  functions  given  by 


X1  x  3 

e  cos  (x2>  +  e  J  sin  x 


1  ’ 


[(xr5)2  +  (x2-4)2  +  (x  -3)23“1/2 


respectively.  The  same  functions  have  been  stu¬ 
died  by  K.  E.  Atkinson  (1980/1982)  using  inte¬ 
gral  equation  methods. 


(iii)  Numerical  Method 


The  numerical  computation  was  done  via  the 
so-called  normal  equations  (using  Cholesky's 
method)  as  described  by  the  author  (cf.  W. 
Freeden  ( 198 3 ) *  chapt.  11). 


The  inner  products  (integrals)  occuring  in 
the  normal  equations  were  computed  by  the  ap¬ 
proximate  sums 


/  f  (  x  ) 
S 


du  a 


1  y 

A  m j  7m2 


1  m. 


) 


where  {x  }  represents  a  (nearly)  uniform 

ml  »m2 

distribution  of  nodes  on  S  with  total  number  A. 

In  our  examples  the  subdivision  (x_  )  on  S 

ml  ,ra2 

is  generated  by  polar  coordinates 


X 

m  |  €  m 

y 

mi  «m2 

^  1 

(sin® 

cosl_  ,  sin® 

sinA  ,cos®  ) 

CM 

E_ 

▼H 

E 

m 

t  m1.m2  m 

1  mi,in2  ml 

with 

ml 

MS 

i  t 

-  m1 

nt  '  2ni 

=  m2 

2s 

t 

CM 

E 

E 

< 

n„  sin(®  ) 

£  m  j 

n_  sin(«  ) 

£ 

E 

V 

o 

i  nl 

,  0  <  m2  <  n2 

^n^,n2:  (fixed)  positive  integersj 

We  choose  especially  nj  =  21  ,  n2  *  42.  Then 
the  total  sum  of  nodes  on  S  is  A  =  562. 

Error  estimates 

We  set 


(1) 

v*  **  -  vlx) 

N 

N 

(2) 

1  at 

v(2)  .  (2) 

We  give  an  impression  of  the  error  for  a  set 
of  selected  points. 


0 

0.4845  E+0 

-0.31  E-6 

-0.34  E-8 

0 

0.1711  E+l 

-0.82  E-4 

-0.79  E-  6 

0 

0.5743  E+0 

-0.64  E-6 

-0.34  E-  S 

0 

0.4388  E+0 

0.11  E-4 

-0.48  E-  7 

0 

0.2422  E+0 

0.84  E-7 

0.12  E-ll 

5 

0.4000  E+0 

0.53  E-8 

0.92  E-  9 

0 

0.2000  E+0 

O.76  E-9 

0.15  E-  9 

Table  2:  Expansion  of  (Surface  2) 


1.0  1.5  2.0 

1.5  0.0  0.0 

3.0  0.0  0.0 

0.0  2.0  0.0 

0.0  4.0  0.0 

0.0  0.0  2.5 

0.0  0.0  5.0 


0.4845  E+O 

-0.91  E-5 

-O.76  E-  6 

0.1711  E+l 

-0.19  E-4 

— 0.72  E—  6 

0.5743  E+0 

-0.13  E-5 

-0.99  E-  7 

0.4388  E+O 

0.54  E-4 

-0.54  E-  6 

0.2422  E+0 

0.35  E-5 

-0.24  E-  7 

0.4000  E+0 

-0.74  E-6 

0.36  E-  8 

0.2000  E+0 

0.19  E-6 

-0.95  E-  9 
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Table 

3: 

Expansion  of 

(Surface  3) 

X1 

X2 

X3 

V(D 

R(l) 

E5 

-(1) 

E6 

0.0 

0.0 

2.0 

0.5000  E+0 

-0.18  E-3 

-0.66  E— 6 

0.0 

0.0 

4.0 

0.2500  E+0 

-0.24  E-3 

0.29  E-6 

0.7 

0.0 

0.0 

0.7375  E+0 

0.82  E-3 

-0.98  E-3 

1.4 

0.0 

0.0 

0.1927  E+l 

-0.35  E-3 

-0.58  E-3 

0.0 

1.0 

0.0 

0.5403  E+0 

0.35  E-2 

-0.52  E-4 

0.0 

2.0 

0.0 

0.4388  E+0 

—0. 12  E-2 

0.17  E-4 

1.0 

1.5 

0.75 

0.7765  E+0 

-0.23  E-2 

-0.39  E-3 

2.0 

3.0 

1.5 

0.3233  E+0 

-0.57  E-3 

-0.91  E-4 

113 


Table  4: 


Expansion  of 


v(2) 


(Surface  1) 


X1 

X2 

X3 

v<2) 

p(2) 

E5 

10 

1.0 

1.5 

2.0 

0.5509  E-l 

0.33  E-  9 

0.23 

E-10 

1.5 

0.0 

0.0 

0.1008  E+0 

0.88  E-  8 

-0.99 

E-  9 

3.0 

0.0 

0.0 

0.4874  E-l 

0.73  E-10 

O' 

a 

0 

1 

E-ll 

0.0 

2.0 

0.0 

0.7352  E-l 

-0.24  E-  8 

-0.20 

E-  9 

0.0 

4.0 

0.3606  E-l 

-0.17  E-10 

-0.10 

E-ll 

0.0 

0.0 

2.5 

0.5788  E-l 

-0.59  E-  9 

-0.35 

E-ll 

0.0 

0.0 

5.0 

0.2862  E-l 

-0.35  E-ll 

-0.51 

E-13 

Table 

5: 

(2) 

Expansion  of  V 

(Surface  2) 

X1 

x2 

X3 

v<2> 

r(2) 

E5 

N 

U 

1.0 

1.5 

2.0 

0.5509  E-l 

0.44  E-  8 

0.12  E-  9 

1.5 

0.0 

0.0 

0.1008  E+0 

-0.44  E-  9 

-0.10  E-  9 

3.0 

0.0 

0.0 

0.4874  E-l 

-0.27  E-  9 

-0.55  E-ll 

0.0 

2.0 

0.7352  E-l 

-0.98  E-  8 

-0.76  E-  9 

0.0 

4.0 

0.3606  E-l 

-0.25  E-  9 

-0.64  E-10 

0.0 

0.0 

2.5 

0.5788  E-l 

0.97  E-  8 

0.13  E-  8 

0.0 

0.0 

5.0 

0.2862  E-l 

O.85  E-  9 

0.98  E-10 
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Appendix  1 : 


maximal  linearly  independent  systems  of  homogeneous 
harmonic  polynomials  for  n  =  3  through  n  =  10  (cf. 
chapter  6). 
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